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The ATLAS Collaboration 



G. Aad^^ B. Abbottiii, J. Abdallahi^, A.A. Abdelalim^^ A. Abdesselami'**, O. Abdinoyi", B. Abi^i^, 

M. Abolins^^ H. Abramowicz'^^ H. Abreu"^ B.S. Acharya"^^''^''^^'', D.L. Adams^^, T.N. Addy^^, 

J. Adelman'^4^ q Adorisio^'^''''^'^'', P. Adragna", T. Adye'^?, S. Aefsky^^, J.A. Aguilar-Saavedrai24b 

M. Aharrouche^i, S.P. Ah[m^\ F. Ahles'*^ A. Ahmad'"^^, H. Ahmed^, M. Ahsan^", G. AieUii33a,i33b^ 

T. Akdogani**'\ T.P.A. Akesson™, G. Akimotoi^^ A.V. Akimov A. Aktas'^^ M.S. Alami, M.A. Alam™, 

S. Albrand^^ M. Aleksa^'', I.N. Aleksandrov^^, C. Alexa^^", G. Alexander'^^^ q Alexandre'^'^, 

T. Alexopoulos'', M. Alhroob^O, M. Aliev'^ G. AUmonti*^^ J. Alison^^o, M. Aliyev'O, P.P AUport^^ 

S.E. AUwood-Spiers^^ J. Almond*^, A. Aloisioi''2'''i02b^ r Alon^™, A. Alonso"^^, M.G. Alviggii02a,i02b^ 

K. Amako^^, C. Amelung^^, A. Amorim'^^''', G. Amoros'^*', N. Amram'^^, C. Anastopoulos'^^, 

T. Andeen^'^, C.F. Anders'^^ K.J. Anderson^", A. Andreazza**''^ V. Andrei5**% X.S. Anduaga™, 

A. Angerami^^, F. Anghinolfi^^, N. Anjos'^"*", A. Annovi'*^, A. Antonaki**, M. Antonelli'*^, 

S. AntonelU'^"''^'', J. Antos^'^'^'', B. Antunovic^^i, R AnulU^^a^ 3 Aoun*^^ G. Arabidze^ I. Aracena''^^ 

Y. Arai^^, A.T.H. Aice^'^, J.P. Archambault^^ S. ArfaouP ", J-F. Arguin''*, T. Argyropoulos^, M. Arikl^^ 

A.J. Armbmster^^, O. Arnaez"*, C. Arnault^'^, A. Artamonov^^, D. Amtinov^'', M. Asai'"*^, S. Asai^^"^, 

R. Asfandiyarov'^^ S. Ask*^, B. Asman^'*^"''"^^'', D. Asner^^, L. Asquith^^, K. Assamagan^^, 

A. Astbury'^^, A. Astvatsatourov^^, G. Atoian^^'*, B. Auerbach^^'*, K. Augsten'^^, M. Aurousseau"*, 

N. Austin", G. Avolio'^^^ r Avramidou^ D. Axen'^^, C. Ay^^^ q Azuelos'^^''', Y. Azuma^^^^, 

M.A. Baak^'', A.M. Bachi"*, H. Bachacoui^*', K. Bachas^'', M. Backes'*'', E. Badescu^S'', P Bagnaia"2a,i32b^ 

Y. Bai^^a^ J Bain'S^ J.T. Baines'^g^ q jg; Baker'^^^ ^ £, Baker^*, S Baker", 

F. Baltasar Dos Santos Pedrosa^^, E. Banas^^ P. Banerjee^^ S. Banerjee'^^ D. Banfi^^'''^^'', A. Bangert^^^, 
V. Bansall^^ S.P. Baranov^"^, S. Baranov^^, A. Barashkou^^, T. Barber^'^, E.L. Barberio^^, 

D. Barberis^O''''^'"', M. Barbero^o, D.Y Bardin*'^ T. Barillari^'^, M. Barisonzi"^ T. Barklow''^^ 
N. Barlow^^, B.M. Barnett'^'', R.M. Barnett'"^, A. Baroncelli'^^a^ A.J. Ban-i'*, F. Barreiro^°, J. Barreiro 
Guimaraes da Costa^'', P. Barrillonii^, R. Bartoldus^'^^ D. Bartsch^", R.L. Bates^^, L. Batkova^'*'*'', 
J.R. Batley^'', A. Battaglia'*', M. Battistin^^, F Bauer'^*', H.S. Bawa'43, M. Bazalova'^S, B. Beare'", 
T. Beau^^ RH. Beauchemin''^ R. Beccherle5"% N. Becerici'**'\ R Bechtle^i, G.A. Beck", H.R Beck'^, 
M. Beckingham'***, K.H. Becks'", A.J. Beddalll**^ A. Beddall^^^ V.A. Bednyakov^^ C. Bee^^ 
M. Begel24, S. Behar Harpaz'^i, RK. Behera^^ M. Beimforde''^, C. Belanger-Champagnei^^ RJ. Bell'^'', 
W.H. Bell"*^, G. Bellai52, L. Bellagambal^^ R Bellina^^, M. Bellomoll9^ A. Belloni^"?, K. Belotskiy^^, 
O. Beltramello^'^, S. Ben Ami'^', O. Benary'^^ Benchekroun"^^ M. Bendel**', B.H. Benedict'''^, 
N. Benekos'*''*, Y Benhammou'^^, G.R Benincasa'^'*'', D.R Benjamin'*'*, M. Benoit**^ J.R. Bensinger^^, 
K. Benslama^^", S. Bentvelsen^*^^, M. Beretta'*^, D. Berge^^, E. Bergeaas Kuutmann'**, N. Berger"*, 
R Berghaus^^^ E. Berglund"*^, J. Beringeri"*, R Bernat^i^, R. Bernhard"*^, C. Bernius'''', T. Berry 
A. Bertini9'».i9b^ y^j Besana^^^-^^'', N. Besson"6, s. Bethke^^, R.M. Bianchi'*^ M. Bianco"'''™, 

0. Biebel^^ J. Biesiada*", M. Biglietti"2aa32b^ jj Bilokon^^, M. Bindii'^" !'^'', S. Binet"^ A. Bingul'^^ 
C. Binii32a,i32b^ Q Biscarat*™, U. Bitenc'^^ K.M. Black^^, R.E. Blair^, J-B Blaiichard"^ G. Blanchot^^, 

C. Blocker22, A. Blondel"*^, W. Blum^i, U. Blumenschein^^, G.J. Bobbink^o^ A. Bocci'*'*, M. Boehler'*!, 
J. Boeki", N. Boelaert^'^, S. Boser", J.A. Bogaerts^^ A. Bogouch'^" *, C. Bohm'^s^ j. Bohmi^s, 

V. Boisvert^^, T. Bold*''^'', V. Boldea^^a, V.G. Bondarenko'''', M. Bondioli'^'^^ Boonekampi^*', 

S. Bordoni", C. Borer*^, A. Borisov*^*, G. Borissov^', I. Borjanovic^^a^ S. Borroni'^^a.nib Bos*"^, 

D. Boscherini*^^, M. Bosman**, H. Boterenbrood*'^^, J. Bouchami^^, J. Boudreau*'^^, 

E. V. Bouhova-Thacker''\ C. Boulahouache^^^, C. Bourdarios*^^, A. Boveia^", J. Boyd^^, I.R. Boyko''^ 

1. Bozovic-Jelisavcic'^'', J. Bracinik*^, A. Braem^^, P. Branchini'^'*^, G.W. Brandenburg^^, A. Brandt', 

G. Brandt^i, O. Brandt^'', U. Bratzler^^S, B. Brau^'^, J.E. Brau""*, H.M. Braun*", B. BreUer*", 
J. Bremer^*^, R. Brenner^^^, S. Bressleri^*, D. Britton^^ R.M. Brochu^'', I. Brock^", R. Brock^*, 
E. Brodet'^^2 q Bromberg^^ G. Brooijmans^"*, W.K. Brooks^"', G. Brown^^^ 

RA. Bruckman de Renstrom^^ D. Bruncko*'*'*'', R. Bruneliere'*^ S. Brunet"**, A. Bruni*^", G. Bruni*'^'', 
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M. Braschil9^ F. Bucci^^, J. Buchanan^i^ P. Buchholzi'*!, A.G. Buckley'*^ I.A. Budagov'^^ B. BudicklO^ 

V. Buscher^i, L. Bugge^'^, O. Bulekov'^'', M. Bunse'*^^ j Buran'^^, H. Burckhart^'', S. Burdin''^ 

T. Burgess'^ S. Burke^29^ £ Busato^^ P. Bussey^^ CP Buszello^''^ F. Butin^^ B. Butler"^-', 

J.M. Butler^i, CM. Buttar^^ J.M. Butterworth^^, T. Byatt^^, J. Caballero^^, S. Cabrera Urban^^^, 

D. Caforioi''^'''"', O. Cakir3^ R Calafiurai"*, G. Calderini^^ R Calfayan^^ R. Calkinsl''^^ L.R Caloba^^^ 

D. Calvet^^ R Camarri"3a,i33b^ j) Cameron"'', S. Campana^'', M. Campanelli", V. Canaleio^a.ioib^ 

F Canelli^°, A. Canepa'5^^ J. Cantero^", L. Capasso'^^aJOib m.D.M. Capeans Garrido^'*, I. Caprmi^^^, 

M. Caprini^s^ M. Capua^^^-^^b, r. Caputo'^?^ ^ Caramarcu^^", R. Cardarem'^Sa^ j CarU^?, 

G. Carlmoi02a, l. Carminati^^^'^^'', B. Caron^ '', S. Caron'*^ G.D. Carrillo Montoya"\ 

S. Carron Montero'^^, A.A. Carter''^, J.R. Carter^^, J. Carvalhoi24a^ Casadei'''^ M.P Casado", 

M. Cascella^^^''''^^'', A.M. Castaneda Hernandez^^^ E. Castaneda-Miranda^^', V. Castillo Gimenez''''', 

N.R Castro'24b, G. Cataldi^^a^ ^ Catinaccio^^, J.R. Catmore^^ A. Cattai^^ G. Cattanii33a,i33b^ 

S. Caughron^^^ j-, (-.^y^i^^a^iesc^ p Cavalleri■'^ D. CavaUi^9^ M. Cavalli-Sforza", V. Cavasimiii22a,i22b^ 

F Ceradmi"'*^'!^'*'', A.S. Cerqueira23=', A. Cerri^^, L. Cerrito''^ F Cerutti'*'', S.A. Cetini**', A. Chafaq'3^^ 

D. Chakraborty"'<'^ K. Chan^, J.D. Chapman^'', J.W. Chapman^^, E. Chareyre''^ D.G. Charlton'^, 

V. Chavda^^ S. Cheatham^i, S. Chekanov^ S.V. Chekulaev'5^^ G.A. Chelkov^^ H. Chen^^, S. Chen32'=, 

X. Chen^'^i, A. Cheplakov^^, V.F. Chepumov^^, R. Cherkaoui El Moursli^^^'^, V. Tchemiatine^^, 

D. Chesneanu^Sa, E. Cheu^, S.L. Cheung^^^, L. Chevalier"^, F Chevallier'^^, V. Chiarella'^^, 

G. Chiefarii''2a,i02b^ £ Chikovani^i, J.T. Childers5^% A. Chilingarov^^ G. Chiodini^^a^ y Chizhov^^ 

G. Choudalakis^", S. Chouridou'^^, I.A. Christidi^^, A. Christov'^^ D. Chromek-Burckhart^'*, M.L. Chu^^", 
J. Chudoba'2^ G. Ciapetti'^^a.iaib^ ^ ^ Ciftci3^ R. Ciftci3^ D. Cinca33, V. Cindro^"*, M.D. Ciobotaru'*'^^ 
C. Ciocca^^^'i^'', A. Ciocioi"*, M. CirilU*'', M. Citterio^^^ A. Clark^^, RJ. Clark'*^ W. Cleland^^^ 

J.C Clemens^^ B. Clement^^ C. Clementi^Sa.wsb^ y. Coadou^^ M. Cobal^^^'' I'^^c^ ^ Coccaro^"^'^"*', 
J. Cochran^'*, J. Coggeshalli'^'*, E. Cogneras'™, A.P Colijn^O^ C Collard"^ N.J. Collins", 

C. Collins-Tooth^^ J. CoUot^^ G. Colon^'^, P. Conde Muino'^^a, E. Coniavitisi^^, M. Consonni^O^, 
S. Constantinescu^Sa, c. Conta"'*^'"'^'', F Conventii'^^a^^ m. Cooke^^, B.D. Cooper"^ 

A. M. Cooper-Sarkar"^ N.J. Cooper-Smith^^, K. Copic^"*, T. CorneHssen^''^'^'"', M. Corradi^'^^ 

F. Corriveau^^ ^ A. Corso-Radu'^^^ ^ Cortes-Gonzalez'^'*, G. Cortiana*^^, G. Costa^^'', M.J. Costa"'^, 

D. Costanzo'3'', T. Costin^", D. Cote''', R. Coura Torres22^ L. Courneyea"^^ G. Cowan^^, C. Cowden^^, 

B. E. Cox^2, K. Cranmer'"^, J. Cranshaw^, M. Cristinziani^", G. Crosetti^^^'^^'', R. Crupi'^^^'''^*', 

S. Crepe-Renaudin^^, C. Cuenca Almenar'^^, T. Cuhadar Donszelmann'^^, M. Curatolo^^, CJ. Curtis'^, 
P Cwetanski*^', Z. Czyczula'^'*, S. D'Auria53, M. D'Onofrio", A. D'Orazio'^'^, C Da Via*^^ 
W. Dabrowski^^ T. Dai*^\ C. Dallapiccola**'*, S.J. Dallison'^'^ *, CH. Daly'^^ M. Dam^^ 

H. O. Danielsson29, D. Dannheim^'', V. Dao'*'^, G. Da^bo50^ G.L. Darlea^^'', W. Davey^^, T. Davidek'^^, 
N. Davidson^^, R. Davidson'^', M. Davies^^ A.R. Davison'^'', I. Dawson"^, R.K. Daya^^, K. De^, 

R. de Asmundis'^^a, s. De Castro''*^ ''^'', PE. De Castro Faria Salgado^", S. De Cecco^^ J. de Graat'^**, 

N. De Groot'^"', P. de Jong"^^ L. De Mora^', M. De OUveira Branco^^, D. De Pedis'^^'', A. De Salvo'-'^^ 

U. De Sanctis'^^^'i^^^ A. De Santo"'^ J.B. De Vivie De Regie"^ G. De Zorzi"2a,i32b^ § Dean'^'^, 

D.V. Dedovich^^ J. Degenhardt'™, M. Dehchar"*^, C. Del Papa"^^^"^^^ J. Del Peso^", 

T. Del Prete'22a''22'', A. Dell'Acqua^^ L. Dell'Asta^^^'^^'', M. Delia Pietra"'2a"', D. della Volpe'^^^'^O^b^ 

M. Delmastro^**, PA. Delsart^^ C. Deluca'"*^, S. Demers'^'*, M. Demichev^^ B. Demirkoz", J. Deng'^^^ 

W. Deng^"^ s p Denisov'^^, J.E. Derkaoui'^^'', F. Derue'^ R Dervan^^, K. Desch^", RO. Deviveiros'^'', 

A. Dewhurst'2^, B. DeWilde^"'', S. DhaUwal'", R. DhulUpudi^^'/, A. Di Ciaccio"'^'''^'', L. Di Ciaccio", 

A. Di Domenico'32'i''32b^ a. Di Girolamo^^, B. Di Girolamo^^ S. Di Luise'^"''''^'"', A. Di Mattia*^^ 

R. Di Nardo'33a.i33b^ a. Di Simone'^^^-'^^'', R. Di Sipio'''^-""', M.A. Diaz^'^ F Diblen'^^ E.B. DiehF, 

J. Dietrich''^ T.A. Dietzsch^*^, S. Diglio"^ K. Dindar Yagci^^, J. Dingfelder^^, C. Dionisi'32a,i32b^ 

R Dita25% S. Dita25% F Dittus^^ F Djama^^ R. Djilkibaev"'^ T. Djobava^', M.A.B. do Vale23% 

A. Do Valle Wemans'^^^ T.K.O. Doan'*, D. Dobos^^, E. Dobson^^ M. Dobson"^^ C Doglioni"^ 

T. Doherty^^ J. Dolejsi'^^, I. Dolenc'^'*, Z. Dolezal'^^, B.A. Dolgoshein^^, T. Dohmae'^", M. Donega^^o, 
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J. Donini^s, J. Dopke"^ A. Doria'O^a^ A. Dos Anjos^i, A. Dottii22a,i22b^ j Dova™, A. Doxiadis^o^ 
A.T. Doyle53, Z. Drasal^^e, M. Dris^ J. Dubbert'^'', E. Duchovnii™, G. Duckeck^^ A. Dudarev^^, 

F. Dudziak'i^ M. Duhrssen L. Duflot'^^ M-A. Dufour**^ M. Dunford^", H. Duran Yildiz^b, 

A. Dushkin22, R. Duxfield^^'^, M. Dwuznik", M. Duren^^^ l Ebenstein'*'*, J. Ebke''^ S. Eckweiler^^ 
K. Edmonds^i, C.A. Edwards^^ K. Egorov", W. Ehrenfeld^^i, T. Ehrich^^, T. Eifert^^, G. Eigeni^, 

K. Einsweileri"*, E. Eisenhandler'^^ T. Ekelof^^^ M. El Kacimi'*, M. EUerti^^, S. EUes'^, F. EUinghaus^i, 
K. Ellis^^, N. EUis^^, J. Elmsheuser^^, M. Elsing^^, D. Emeliyanov'^^, R. Engelmann'"*^, A. Engl^**, 

B. Epp^^, A. Eppig^^, J. Erdmann^'^, A. Ereditato'^, D. Eriksson i'^^", I. Ermoline^^ J. Emst\ M. Ernst^"*, 
J. Emwein"6, d. Erredei^, S. Erredei^'*, E. Ertel^i, M. EscaUer"^ C. Escobar^'^^^ x. Espinal Curallii, 

B. Esposito^^, A.I. Etienvre'^^, E. Etzion^^^ H. Evans", L. Fabbrii^"''^'', C. Fabre^'', K. Facius^^ 

R.M. Fakhmtdinovi28, S. Falciano'^^a^ Y. Fang'^i, M. Fanti**'^" **'^'', A. Farbin^ A. Farillai34a^ j. Farley^^?, 
T. Farooque'^^, S.M. Farrington^'^, R Farthouat^^, R Fassnacht^^, D. Fassouliotis^, B. Fatholahzadeh^^^, 
L. Fayard'15, F. Fayette^'*, R. Febbraro^^ P. Federic''*'*^ O.L. Fedini^i, w. Fedorko^^, L. Feligioni^^ 

C. U. Felzmann^^ C. Feng^^d^ g j pgngSO^ ^ g Fenyuk^^s^ j Ferencei^'*^'', J. Ferland^^ B. Femandes^^^a, 
W. Fernando'"^, S. Ferrag^^, J. Ferrando^'^, V. Ferrara^', A. Ferrari"'^, R Ferrari'"^, R. Ferrari"^^ 

A. Ferrer"''', M.L. Ferrer'*^, D. Ferrere'*^, C. Ferretti^^, M. Fiascaris^'^, F. Fiedler^^ A. Filipcic^"*, 

A. Filippas^ F Filthautio^, M. Fincke-Keeleri^^ M.C.N. Fiolhaisi^^a, l. FioriniH, A. Firan^^, 

G. Fischer^i, M.J. Fisher'"'^, M. Flechli'^^ I. Fleck^^i, J. Fleckner^i, R Fleischmann"^^ § Fleischmann^", 
T. Flick"^ L.R. Flores Castilloi^\ M.J. Flowerdew'^'', T. Fonseca Martin^'', A. Formica'^*', A. Forti^^^ 

D. Fortin'^*", D. Fournier^'^, A.J. Fowler'^'*, K. Fowler'^'', H. Fox^\ R Francavilla'22aa22b^ 

S. Franchino'i'^'''"^'', D. Francis^^, M. FrankUn^'', S. Franz^^, M. FraternaU'!'^'''!'^'', S. Fratina^^O, 
J. Freestone*^, S.T. French-^^, R. Froeschl^^, D. Froidevaux^^, J.A. Frost^^, C. Fukunaga^^^, 

E. Fullana Torregrosa^, J. Fuster^^*', C. Gabaldon^", O. Gabizon'™, T. Gadfort^^, S. Gadomski"*^, 

G. Gagliardi50'''5'^^ P. Gagnon", C. Galea''^ E.J. Gallas'^^ V. Gallo'^ B.J. Gallopi^^, R Gallus'^s, 

E. Galyaev'*'^, K.K. Gan^^^ Y.S. Gao^'*^-^ A. Gaponenko^^, M. Garcia-Sciveresi'*, C. Garciai^^ J.E. Garcia 
Navarro49, R.W. Gardner^", N. Garelli^'', H. Garitaonandia'^^ V. Garonne^^, C. Gatti^^, G. Gaudioll'^^ 
V. Gautard'^*', R Gauzzii32a,i32b^ j l Gavrilenko'^^^ q Gay^^\ G. Gaycken^", E.N. Gazis^ R Ge32d, 

C. N.R Geei29, Ch. Geich-GimbeP, K. Gellerstedti'^^^ i'^Sb^ c. Gemme5°^ M.H. Genest*^**, 

S. Gentile"2a,i32b^ p Qeorgatos^ S. George^*', A. Gershoni52, h. Ghazlane^^Sd, n. Ghodbane^^, 

B. Giacobbel9^ S. Giagu"2a.i32b^ y. Giakoumopoulou^ V. Giangiobbei22a,i22b^ p Gianotti29, 

B. Gibbard24, A. Gibson^^?, S.M. Gibson''**, L.M. Gilbert''^ M. Gilchriesei", V. Gilewsky*^', 

D. M. Gingrich2 '', J. Ginzburgi52^ n. Giokaris**, M.R Giordani"^^'' i'^^", R. Giordano' "2a, i02b^ pM. Giorgi'^ 
R Giovannini'^'', RF Giraud2'', R Girtler'^2^ £, Giugni**'^'', R Giustii'''\ B.K. Gjelsten"^ L.K. Gladilin*^^, 

C. Glasman^o, A. Glazov'", K.W. GUtza"^ G.L. Glonti^^ J. Godfrey'''2, j. Godlewski29, M. Goebel^', 
T. Gopfert*^ C. Goeringer^', C. G6ssling42, T. Gottfert^^ V. Goggi"^^''!^''''', S. Goldfarb^'^, D. Goldin^^, 
T. Golling'^'', A. Gomes'24a^ L.S. Gomez Fajardo"", R. Gongalo^^, L. Gonella2'', C. Gong^2b^ 3 Gonzalez 
de laHoz"'^, M.L. Gonzalez Silva26, S. Gonzalez-Sevilla'*^ J.J. Goodson'"*^, L. Goossens2^ 

H. A. Gordon24, 1. Gorelov'°^ G. Gorfine^'^^ B. Gorini29, E. Gorini^2a,72b^ ^ Gorisek'^'', E. Gomicki^^ 
B. Gosdzik"", M. Gosselink'"^, M.l. Gostkin^^, 1. Gough Eschrich"'2, M. Gouighri"^^ D. Goujda^li'2^^ 
M.R Goulette^'', A.G. Goussiou'^**, C. Goy^, 1. Grabowska-Bold"'2'', R Grafstrom2'', K-J. Grahn"^^, 

S. Grancagnolo'^ V. Grassi"^^ V. Gratchev'2i, N. Grau^"*^ Gray^"*,;^ j^^^ Gray'''^ E. Graziam'3'*^ 

B. Green'^, T. Greenshaw^^ Z.D. Greenwood24'/, I.M. Gregor"", P. Grenier"'^, E. Griesmayer'*^, 
J. Griffiths'3^ N. Grigalashvili^^ A.A. Grillo'^"^, K. Grimmi'*^, S. Grinstein", Y.V. Grishkevich^^, 

M. Groh'^'^, M. Groll^', E. Gross'™, J. Grosse-Knetter^'^, J. Groth-Jensen^^, K. Grybel"", C. Guicheney^^, 
A. Guida^2a,72b^ j OuiUemin'*, H. Guler^^ j Gunther'25, B. Guo'", A. Gupta^^, Y. Gusakov^^ 
A. Gutierrez^^ R Gutierrez'", N. Guttmani52^ q Gutzwiller"\ C. Guyot'^*, C. Gwenlan'^^, 

C. B. Gwilliam", A. Haas'^^, S. Haas2^ C. Haber'", H.K. Hadavand^^ D.R. Hadley'^, R Haefner'^^ 

R. Hartel'^^ Z. Hajduk^^, H. Hakobyan'^^ J. Haller'" ^ K. Hamacher'^^ A. Hamilton^'', S. Hamilton'^o^ 
L. Haii32i', K. Hanagaki'16, M. Haiice'20, C. Handel^\ P. Hallke5^^ J.R. Hansen^^, J.B. Hansen^^, 
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J.D. Hansen^^ P.H. Hansen^^, T. Hansl-Kozanecka""^, P. Hanssoni^S, K. Hara^^S, Q.A. Hare""^, 

T. Harenberg^^^, R.D. Harrington^', O.M. Harris'^^, K Harrison'"', J. Hartert"'^, F. Hartjes'"^, A. Harvey^^, 

S. Hasegawa"*', Y. Hasegawa"'", K. Hashemi^^, S. Hassani'^*', S. Haug'*', M. Hauschild^'', R. Hauser^*, 

M. Havranek'25, c.M. Hawkes'^, R.J. Hawkings^^, T. Hayakawa*'^, H.S. Hayward^^ S.J. Haywood'^^, 

S.J. Head^2^ y. Hedberg™, L. Heelan^^, S. Heim^^ B. Heinemann'^, s. Heisterkamp^^, L. Helary'', 

M. Heller"^ S. Hellmani^Sa.wsb^ q Helsens", T. Hemperek^o, R.C.W. Henderson''^, M. Henke^*^ 

A. Henriciis^^, A.M. Henriques Correia^^, S. Henrot-Versille"^, C. Hensel^^, T. Henfi'^^, Y. Hernandez 

Jimenez""^', A.D. Hershenhorn'^', G. Herten"***, R. Hertenberger^^, L. Hervas^^, N.R Hessey'*^^, 

E. Hig6n-Rodriguezi66, j.c. HiU^'', K.H. Hiller^i, S. Hillerti''5'''i''5*', S.J. Hillier''', I. Hinchliffe^'', 

E. Hines'20, M. Hirose'"^, F. Hirsch42, D. Hirschbuehl"^ J. Hobbs'''^ N. Hod'^z, M.C. Hodgkinson'^^, 

R Hodgson'39, a. Hoecker^^, M.R. Hoeferkamp'^^ J. Hoffman^^, D. Hoffmann^^ M. Hohlfeld*^', 

T. Holy'27, J.L. Holzbauer^*^, Y. Homma*'\ T. Horazdovsky'^^, T. Hori'^^, C. Horn'43, s. Horner^^, 

S. Horvat^^, J-Y. Hostachy^^, S. Hou'^O, A. Hoummada'^Sa, t. Howe^'', J. Hrivnac"^ T. Hryn'ova'', 

RJ. Hsu''^'', S.-C. Hsu''', G.S. Huang'", Z. Hubacek^^?, p. Hubaut^^ R Huegging^", E.W. Hughes3^ 

G. Hughes^', M. Hurwitz^", U. Husemann"", N. Huseynov'", J. Huston^^, J. Huth^^, G. Iacobucci'°2^ 

G. lakovidis^, I. Ibragimov'^', L. Iconomidou-Fayard"^, J. Idarraga'^^'', P. lengo'*, O. Igonkina"^^, 
Y. Ikegami^^, M. Ikeno^^, Y. nchenko^^ D. Iliadisl5^ T. Incel^^ R loannou^ M. Iodice'3'*^ 

A. Irles Quiles'^*', A. Ishikawa^'', M. Ishino''^, R. Ishmukhametov^'*, T. Isobe'^", V. Issakov'^"*'*, 

C. Issever"^ S. Istin'^% Y. Itoh"", A.V. Ivashin'^^, W. Iwanski3^ H. Iwasaki'^^ J.M. Izen'^", V. Izzo'^^a^ 

B. Jackson'20, J.N. Jackson^^ R Jackson'^^, M.R. Jaekel^'', V. Jain", K. Jakobs"*^, S. Jakobsen^^, 
J. Jakubek'27, D.K. Jana'", E. Jansen'04, a. Jantsch^^, M. Janus''^ R.C. Jared'^', G. Jarlskog™, 
L. Jeanty^"^, I. Jen-La Plante^", R Jenni^^, R Jez^^, S. Jezequel'', W. Ji''^ J. Jia''*'', Y. Jiang^^b^ 

M. Jimenez Belenguer^^ S. Jin^^a^ q jinnouchi'^^ D. Joffe^^^ Johansen'^'^^ '^'Sb, K.E. Johansson'^^^, 
R Johansson'3'^, S Johnerf*', K.A. Johns^ K. Jon-And'^^^'^''^*', G. Jones*^^ R.W.L. Jones^', T.J. Jones", 
RM. Jorge'24a, j. Joseph''^, V. Juranek'^^, R Jussel^^, V.V. Kabaciienko'^^, M. Kaci'^^, A. Kaczmarska'^ 
M. Kado"^ H. Kagan'"'', M. Kagan^^ S. Kaiser^^, E. Kajomovitz'^i, S. Kalinin'", L.V. Kalinovskaya^^ 
A. Kalinowski'^", S. Kama^', N. Kanaya'^'', M. Kaneda'^^, v.A. Kantserov^*', J. Kanzaki'^^ B. Kaplan'^"*, 
A. Kapliy^", J. Kaplon^^, D. Kar'*^ M. Kara gounis^**, M. Karagoz Unel"**, V. Kartvelishvili^', 
A.N. Karyukhin'28, L. Kashif", A. Kasmi^^, R.D. Kass'"'', A. Kastanas'^, M. Kastoryano'^'', M. Kataoka", 
Y. Kataoka'^'', E. Katsoufis^, J. Katzy^i, V. Kaushik^, K. Kawagoe^"^, T. Kawamoto'^'', G. Kawamura*', 
M.S. Kayl'O^ F Kayumov'^'', V.A. Kazanin'^^, M.Y Kazarinov^^ J.R. Keates^^ R. Keeler"^^ 
RT. Keener'20, R. Kehoe^*^, M. Keil^", G.D. Kekelidze^^ M. Kelly**^ M. Kenyon^^ O. Kepka'^^, 
N. Kerschen^'', B.R Kersevan^^, S. Kersten'", K. Kessoku'^'', M. Khakzad^**, F. Khalil-zada"\ 

H. Khandanyan'64^ ^ Khanov"^, D. Kharchenko''^ A. Khodinov'^'', A. Khomich^^a, G. KhoriauU^o, 
N. Khovanskiy65 V. Khovanskiy'*^ E. Khramov^^ J. Khubua^^, H. Kim'', M.S. Kim^, RC. Kiffli^S, 

S.H. Kim™, O. Kind'5, R Kind'", B.T. King", J. Kirk'^^, G.R Kirsch"*^, L.E. Kirsch^^, A.E. Kiryunin^^, 

D. Kisielewska^^ T. Kittelmann'", H. Kiyamura^^ E. Kladiva''*"'', M. Klein", U. Klein", 

K. Kleinknecht^i, M. Klemetti^^ A. Klier'™, A. Klimentov^", R. Klingenberg'*^, E.B. Klinkby'^"^ 

T. Klioutchnikova^'^, RF Klok'"", S. Klous"'^ E.-E. Kluge5**^ T. Kluge", R Kluit"'^ M. Klute^"^ 

S. Kluth'*'^, N.S. Knecht'5\ E. Kneringer^^ B.R. Ko^^, T. Kobayashi'^^, M. Kobel''^ B. Koblitz^^ 

M. Kocian''*^ A. Kocnar"^ R Kodys'^^, K. Koneke^', A.C. Konig"'^^ ^ Koenig*^', L. Kopke*^', 

R Koetsveld"''', R Koevesarki^o, T. Koffas^", E. Koffeman'^s, R Kohn^"^ 2. Kohout'^^, T. Kohriki^^, 

H. Kolanoskii^ V. Kolesnikov^^ I. Koletsou", J. Koll*^ D. KoUar^^, S. Kolosi^^'', S.D. Kolya*^ 

A.A. Komar*^"*, J.R. Komaragiri''*^ T. Kondo^*^, T. Kono"*' ^ R. Konoplich"'^ S.R Konovalov^'', 

N. Konstantinidis^^, S. Koperny", K. Korcyl^^ K. Kordas'^^ A. Korn'", I. Korolkov", E.V. Korolkova'^^, 

V.A. Korotkov'2^ O. Kortner''^ R Kostka-*', V.V. Kostyukhin^o, S. Kotov^^, VM. Kotov*'^ K.Y Kotov'"'', 

C. Kourkoumelis^, A. Koutsman'"^, R. Kowalewski'^^, H. Kowalski^', T.Z. Kowalski^^, W. Kozanecki'^^, 
A.S. Kozhin'28, V. Kral'^'', V.A. Kramarenko^^, G. Kramberger^'^, M.W. Krasny^^, A. Krasznahorkay"^*, 
A. Kreisel'52^ p Krejcii^'', J. Kretzschmar'^^ N. Krieger^"^ p Krieger^^'', K. Kroeninger^"^ jj Kroha^^, 
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J. Kroll'2'\ J. Kroseberg2", J. Krstic'^a, u. Kruchonak6^ H. Krtiger^o, Z.V. Krumshteyn*'^ T. Kubota^^'*, 
S. Kuehn^^ A. Kugel^**^ T. Kuhl'", D. Kuhn*'^, V. Kukhtin'^^ Y. Kulchitsky''", S. Kuleshov^"', 
C. Kummer''*^, M. Kuna*^^, J. Kunkle'^", A. Kupco^^s^ j| Kurashige''^, M. Kurata'^'^, L.L. Kurchaninov'^**^ 
Y.A. Kurochkin''", V. Kus'^s, r. Kweel^ L. La Rotonda^*"'''^'^'', J. Labbe"*, C. Lacasta'f'^ F. Lacavai32a,i32b^ 
H. Lacker^^, D. Lacour^^, V.R. Lacuesta^*''', E. Ladygin*'^, R. Lafaye"*, B. Laforge^*, T. Lagouri^O, S. Lai'*^, 
M. Lamanna^^, C.L. Lampen^ W. Lampl*', E. Lancon"6, u. Landgraf***, M.RJ. Landon^^^ j l Lane*^, 
A.J. Lankford'''^, F. Lanni^^, K. Lantzsch^'', A. Lanza^''^^ S. Laplace^, C. Lapoire^^ J.E Laporte^^^ 
T. Lari^^", A. Earner' M. Lassnig^^, P. Laurelli'^^, W. Lavryseni"*, P. Laycock", A.B. Lazarev^^, 
A. Lazzaro*^^'*^'', O. Le Dortz''^ E. Le Guirriec^^ E. Le Menedeu'^^, M. Le Vine^"*, A. Lebedev^'*, 

C. Lebel'^^ T. LeCompte^ R Ledroit-Guillon^^, H. Lee^o^ J.S.H. Lee'^^, s.C. Leei^", M. Lefebvrei^^ 
M. Legendre'^'', B.C. LeGeyt'^o^ p Legger^**, C. Leggett'"*, M. Lehmacher^'', G. Lehmann Miotto^^, 
X. Lei'', R. Leitner'26, D. Lellouch'^", J. Lellouch^^ V. Lendermann5**^ K.J.C. Leney^^ T. Lenz^'^^ 
G. Lenzen"^, B. Lenzi'^*', K. Leonhardt'*^ C. Leroy'^^ J-R. Lessa^d'^^ C.G. Lester^^ 

A. Leung Fook Cheong"\ J. Leveque*^ D. Levin*^, L.J. Levinsoni™, M. Leyton''^, H. Li"\ S. Li^K 
X. Z. Liang3^ Z. Liangi^o ™, B. Libertii33a^ p. Lichard^^, M. Lichtnecker''^ K. Lie'*"*, W. Liebigi^^ 
J.N. Lilley", H. Lim^ A. Limosani'^*', M. Limper*'^ S.C. Lin'^°, J.T. Linnemann^^ E. Lipeles'^O^ 
L. Lipinskyi25, a. Lipniackai^ T.M. Liss'^"*, D. Lissauer^^, A. Lister"*^, A.M. Litke"^, C. Liu^^, 

D. Liu^^^ ", H. Liu**^ J.B. Lm^\ M. Liu^^b^ j pj^sg^ y. Liu^^b^ jyj pivan'!'''' !''"', A. Lleres^^ 
S.L. Lloyd^^ E. Lobodzinska^i, R Loch*', W.S. Lockman'^^, S. Lockwitz'^"^, T. Loddenkoetter^", 

RK. Loebinger**2^ ^ Loginov'^'*, C.W. Loh'*'^, T. Lohse'^ K. Lohwasser'*^ M. Lokajicek'^s, R.E. Long^^ 
L. Lopesi24a^ D. Lopez Mateos^'^'', M. Losada'*'!, R Loscutoffi"*, X. Lou'^o, A. Lounis'^^ K.R Loureiro^o^ 
L. Lovasl'^^ J. Love^i, RA. Love^\ A.J. Lowe^\ R Lu^^, H.J. Lubattii^^ C. Luci"2a,i32b ^ Lucotte^^, 
A. Ludwig43, D. Ludwig^i, L Ludwig^^ R Luehring", L. Luisa^''^''''^^^'^ , D. Lumb^^ L. Luminari"2a^ 

E. Lund^'^, B. Lund-Jensen'^^, B. Lundberg^^, J. Lundberg^^, J. Lundquist^^, D. Lynn^^, J. Lys'^, 

E. Lytken™, H. Ma^^, L.L. Ma"\ J.A. Macana Goia^^ G. Maccarrone'*'', A. Macchiolo'^^ B. Macek'''*, 

J. MachadoMiguens'^"*^, R. Mackeprang^^, R.J. Madaras'"*, W.R Mader"*^, R. Maenner^^*^, T. Maeno^"^, 

R Mattig'^^ S. Mattig'*', RJ. Magalhaes Martins'^^a^ p. Magradze^', Y. Mahalalel'^^, K. Mahboubi^^, 

A. Mahmoodi, C. Maiani'^^ajsib^ q Maidantchik^^", A. Maio'24a^ s. Majewski^"*, Y. Makida'^*', 

M. Makouskii28, N. Makovec^'^ Ra. Malecki^^, R Malecki^^, V.R Maleev'^i, R Malek^^, U. Mallik^^^ 

D. Malon^, S. Maltezos^, V. Malyshev^"'^, S. Malyukov^^, M. Mambelli^", R. Mameghani^*, J. Mamuzic'*\ 

L. Mandelli^^", L Mandic''"^, R. Mandrysch'^, J. Maneira'2'^^ RS. Mangeard^**, ID. Manjavidze''^, 

P.M. Manning'^^, A. Manousakis-Katsikakis**, B. Mansoulie'^*', A. Mapelli^', L. Mapelli^'', L. March 

J.F. Marchand'*, F. Marchese'^'^''- G. Marchiori^^ M. Marcisovsky'^^, CP. Marino''', F. Marroquim2^^ 

Z. Marshall^^,;^ § Marti-Garcia'^'', A.J. Martin", A.J. Martin'"''*, B. Martin^^ B. Martin^^ RR Martin'^O, 

J.R Martin^^, T.A. Martin''^, B. Martin dit Latour"*^, M. Martinez", V. Martinez Outschoom^'^, 

A. Martini"''^, A.C. Martyniuk**^, F. Marzano'^2^ A. Marzin'^^, L. Masetti^", T. Mashimo'^'', 

R. Mashinistov^'', J. Masik**^, A.L. Maslennikov'"'', L Massa'^"''^'', N. Massol"', A. Mastroberardino^'''''^''^ 

T. Masubuchi'^'', P Matricon"^, H. Matsunaga'^"*, T. Matsushita^'^, C. Math-avers"^ ", S.J. Maxfield^^ 

A. Mayne'^^ R. Mazini'^", M. Mazur''^ M. Mazzanti**'^^ J. Mc Donald^^ S.P Mc Kee^^, A. McCam'^'*, 

R.L. McCarthy'''^ N.A. McCubbin'^'', K.W. McFarlane^'', H. McGlone^^ G. Mchedlidze^', 

S.J. McMahon'29, R.A. McPherson"^**'^ A. Meade^'', J. Mechnich'^^, M. Mechtel'", M. Mediums'*', 

R. Meera-Lebbai'", T.M. Meguro"^, S. Mehlhase'", A. Mehta^^ K. Meier5*^^ B. Meirose'*^, 

C. Melachrinos^", B.R. Mellado Garcia''^', L. MendozaNavas'^', Z. Meng'^^" '\ S. Menke^^, E. Meoni", 
R Mermod"^ L. Merola"'2aao2b^ q Meroni**'^^ RS. Merritt^", A.M. Messina^^ J. Metcalfe"'^ 

A.S. Mete*^"^ j.p Meyer'^^, J. Meyer'^^^ j Meyer^^, T.C. Meyer^^, W.T. Meyer^"^ j Miao^^d^ g Michal^^, 
L. Micu25^ R.R Middleton'2^, S. Migas'^^, L. Mijovic'^'*, G. Mikenberg'""^, M. Mikestikova'^^, 
M. Mikuz^"^ p,;^ Miller'43, W.J. Mills CM. Mills^^, A. Milov'™, D.A. Milstead'''5a.i45b^ 

D. Milstein'™, A.A. Minaenko'^^, M. Minano'^^, LA. Minashvili''^, A.I. Mincer'"**, B. Mindur^^, 
M. Mineev^^ Y. Ming'^"^ p ]y[ ^^n^ q Mirabelli'^^a, s. Misawa^"*, S. Miscetti'*'', A. Misiejuk'^^, 
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J. Mitrevski"?, V.A. Mitsou'^^, P.S. Miyagawa*^^ J.U. Mjornmark''^ D. Mladenov^^, T. Moa^"^^^'^*^^, 
S. Moed", V. Moeller^^ K. Monig^i, N. Moser^", W. Mohr^**, S. Mohrdieck-Mock'^^ R. Moles-VaUsi^^, 
J. Molina-Perez^^, J. Monk^^, E. Monnier^^, S. Monlesano**'^'^ **'^'', F. Monticelli^'', R.W. Moore^, 

C. Mora Herrera'^^, A. Moraes^^, A. Morais'^"*", J. Morel^'*, G. Morello^'''''^'''', D. Moreno^''^ M. Moreno 
Llaceri^*^, P. Morettini^'''', M. Morii^^, A.K. Morley^^, G. Mornacchi^'', S.V. Morozov^*^, J.D. Morris^^, 
H.G. Moser^^, M. Mosidze'i, J. Moss^''^, R. Mounti'^^ E. Mountricha"^, S.V. Mouraviev'^'*, 

E. J.W. Moyse**-*, M. Mudrinici^b, R Mueller5**^ J. Mueller'23, K. Mueller^", T.A. Muller'^^ 

D. Muenstermann'*^, A. Muir"'^, Y. Munwes^^^, R. Murillo Garcia'^^, W.J. Murray'^^, I. Mussche'^^, 

E. Mustoi02a,i02b^ ^ q Myagkov'^s, M. Myska'^s, j. Nadal", K. Nagai^^^, K. Nagano^*', Y. Nagasaka^o, 
A.M. Nairz29, K. Nakamurai^"*, I. Nakano'i"*, H. Nakatsuka^^ G. Nanava^o, A. Napier^'^'', M. Nash"'"?, 
N.R. National, T. Nattermann^o, T. Naumann^i, G. Navarroi^i, S.K. Nderitu^o, H.A. Neal^^, E. Nebot^*^, 
P Nechaeva''^, A. Negri' ^'^^''^i^'', G. Negri^^ A. Nelson'^'^, T.K. Nelson'^^, s. Nemecek'^S, R NemethylO^ 
A.A. Nepomuceno^^^, M. Nessi^^, M.S. Neubaueri^'*, A. Neusiedl^\ R.M. Neves'^^ P Nevski^'*, 

F. M. Newcomer'20^ g Nickerson''^, R. Nicolaidou'^^, L. Nicolas''^, G. Nicoletti'*'', B. Nicquevert^^, 

F. Niedercorn"^, J. Nielsen'^'', A. Nikiforov'^, K. Nikolaev^^, I. Nikolic-Audit''**, K. Nikolopoulos^, 

H. Nilsen"*^, P. Nilsson^, A. Nisati'^^'', T. Nishiyama^^, R. Nisius^^, L. Nodulman^, M. Nomachi'^^, 

I. Nomidisi^^ M. Nordberg^'', B. Nordkvisti''^'''i''5'', D. Notz^i, J. Novakova'^e, M. Nozaki^^, 
M. Nozicka"*', I.M. Nugent'^^^ A.-E. Nuncio-Quiroz^'', G. Nunes Hanninger^", T. Nunnemann'^, 
E. Nurse^\ D.C. O'Neil'42, v. 0'Shea5^ EG. Oakham^*^ H. Oberlack'''', A. Ochi'^^ S. Oda'^^, 

S. Odaka^^, J. Odier^^ H. Ogren", A. Oh^^ S.H. Oh'*"^ q q Ohm'^Sa.wsb^ j Ohshimai'^i, H. Ohshita''^'', 
T. Ohsugi5^ S. Okada^'', H. Okawa''^^ Y. Okumurai'^\ T. Okuyama'^'^, A.G. 01chevski''5, M. Oliveirai^^a, 
D. Oliveira Damazio^, J. Oliver^'', E. Oliver Garcia'^^, D. Olivitoi^", A. Olszewski^^, J. Olszowska^^, 

C. Omachi'^^ ^ A. Onofre'^^a^ PU.E. Onyisi^", C.J. Oraml5^^ M.J. Oreglia^", Y. Oreni52, 

D. Orestano^^'^^ i^^b^ j Orlov^^^ C. Oropeza Barrera^^ R.S. Orr^^^ E.G. Ortega^^^, B. Osculati^^'''^'^'', 

R. Ospanov'20, c. Osuna'\ J.P Ottersbachi°^ R Ould-Saada^i'', A. Ouraou^^^, Q. Ouyang'^a^ Owen*^^ 

S. Owen'^*^, A Oyarzun^"', V.E. Ozcan'''', K. Ozone^^, N. Ozturk'', A. Pacheco Pagesi\ 

C. Padilla Aranda'i, E. Paganisi^*^, C. Pahl'^^ R Paige^^, K. Pajchel^i^ S. Palestini^^, D. Pallin^^ 

A. Palma^24a^ j Palmer^^, YB. Pan"\ E. Panagiotopoulou*^, B. Panes^''', N. Panikashvili**^, 

S. Panitkin24, D. Pantea2^^ M. Panuskova^^s, v. Paolone'^^ Th.D. Papadopoulou'', S.J. Park^^^ 

W. Park24.^, M.A. Parker^', S.I. Parker^^, R Parodi^''^'^''*', J.A. Parsons^'^, U. Parzefall'*^ 

E. Pasqualucci'^^a^ A. Passerii^'^% R Pastore'^'^^'i^^b^ Pastore^^, G. Pasztor S. Pataraia^^, 
J.R. Pater**^ S. Patriceni"'2a.i02b^ ^ Patwa^^, T. Pauly^^, L.S. Peak'^^, m. Pecsy'^^^a^ 

M.I. Pedraza Morales'^' , S.V. Peleganchuk"", H. Peng'^^ A. Penson^"*, J. Penwell^^ M. Perantoni^^^, 

K. Perez^"* ', E. Perez Codina'\ M.T. Perez Garcfa-Estafli*^'', V. Perez Reale^^^ l Perini^''^'^'"', 

H. Pernegger^^, R. Perrino''^*, S. Persembe^^, P. Perus^'^, V.D. Peshekhonov^^, B.A. Petersen^^, 

T.C. Petersen^^ E. Petit**^ C. Petridou'^^ E. Petrolo'^^a^ p Petrucci'^'*^ '^'*'', D Petschull'*^ M. Petleni^'*^^ 

R. Pezoa^'*', A. Phan'^^ A.W. PhiUips^^ G. Piacquadio^'', M. Piccinini!'^'''!'^'', R. Piegaia^^, J.E. Pilcher^o, 

A.D. Pilkington*2^ j pjna^a^ m. Pinamontii^2'''i^3c j l Pinfold^, B. Pinto™^ C. Pizio^^^'^^*', 

R. Placakyte"^', M. Plamondon'*'^ M.-A. Pleier^'^, A. Poblaguev'^'*, S. Poddar^^^ R Podlyski^^ 

P Poffenbergeri*'*^, L. Poggioli'i^ m. Pohl'*'^, R Polci^^ G. Polesellol''^^ A. Policicchio'^**, A. Polinii^^, 

J. PoU^^, V. Polychronakos^'*, D. Pomeroy^^, K. Pommes^^, P. Ponsot'^*', L. Pontecorvo'^^a^ g q Pope*^, 

G. A. Popeneciu25^ D.S. Popovic'2^ A. Poppleton^", J. Popule^^^, x. Portell Bueso'^^ R. Porter^^^^ 
G.E. Pospelov^^, S. Pospisil'^?, M. Potekhin^^, I.N. Potrap''^ C.J. Potterl'^^ C.T. Potter*^ K.P Potter^^^ 
G. Poulard^^, J. Poveda'^', R. Prabhu^", R Pralavorio^^ S. Prasad^^, R. Pravahan^, L. PribyP, D. Price^\ 
L.E. Price^, P.M. Prichard^^ D. Prieur'^^ M. Primavera^^a^ Prokofiev^^, R. Prokoshin^^'', 

S. Protopopescu^^, J. Proudfoot^, X. Prudenf*^, H. Przysiezniak"*, S. Psoroulas^", E. Ptacek^'"*, 

C. Puigdengoles'^ J. Purdham^'', M. Purohit^^ P Puzo''^ Y. Pylypchenko'^^, M. Qi^^c^ j Qi^&l^ 

W. Qiani29^ ^ q^^41^ ^ Quadt^'^, D.R. Quarrie''^, W.B. Quayle^^^ R Quinonez3l^ M. Raas^**, 

V. Radeka24, V. Radescu^^'', B. Radics^^, T. Radorl«^ R Ragusa^^"-*^'', G. Rahal"^ A.M. Rahimiio^, 
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S. Rajagopalan'^'*, M. Rammensee^**, M. Rammes^^\ F. Rauscher'**, E. Rauter^^, M. Raymond'^^, 

A.L. Readily D.M. Rebuzzi'!''''''!!'^'', A. Redelbachi^^ G. Redlinger^^, R. Reece'^o, K. Reeves^", 

E. Reinherz-Aronis'^2^ ^ Reinsch'^'^, I. Reisinger^^, D. Reljic^^a^ q Rembser^'', Z.L. Ren^^'', R RenkeP^, 

S. Rescia^"*, M. Rescigno'^^a^ 5 Resconi^^^, B. Resende'^*", R Reznicek^^e^ ^ Rezvani'", A. Richards'^'', 

R.A. Richards^^ R. Richter^^, E. Richter-Was^^-", M. Ridel^^ M. Rijpstra'°5, M. Rijssenbeeki'*'', 

A. Rimoldiii^'*'"^'', L. Rinaldil^^ R.R. Rios^^, I. Riu", F. Rizatdinova"^, E. Rizvi^^ 

D.A. Roa Romero'^', S.H. Robertson*^'' A. Robichaud-Veronneau'*^ D. Robinson^^, JEM Robinson^^, 

M. Robinson! I'*, A. Robson^^ J.G. Rocha de Lima^^^^, C. Rodai22a,i22b^ r^^j^ j^^^ Santos^^, 

D. Rodriguez'^\ Y. Rodriguez Garcia^^, S. Roe^^, O. R0hneii'^, V. Rojo\ S. Rollii^°, A. Romaniouk^^, 

V.M. Romanov''^ G. Romeo^*', D. Romero Maltrana^''', L. Roos^^ E. Ros'^*', S. Rosati"^ 

G. A. Rosenbaumi", L. Rosselet'*'', V. Rossetti^^ L.R Rossi50^ M. Rotaru^^a, j. Rothberg'^^ 

D. Rousseau"^ C.R. Royon'^*', A. Rozanov^^ Y. Rozen'^', X. Ruan'^^ B. Ruckert'^^ N. Ruckstuhli''^, 
V.I. Rud^'', G. Rudolph^^^ p RyhrSSa^ p. Ruggierii^"*", A. Ruiz-Martinez^"^, L. Rumyantsev*'^ 
Z. Rurikova"^*, N.A. Rusakovich*^, J.R Rutherfoord^, C. Ruwiedel^", R Ruzicka^^^, YF. Ryabov'2\ 
R Ryan**^ G. Rybkin"^ S. Rzaeva'", A.F Saavedra"*'^, H.F-W. Sadrozinski"^, R. Sadykov''^ 

H. Sakamoto'^'*, G. Salamanna"'^ A. Salamon'^Sa^ jyj 3 Saleem'^i, D. SaUhagic^^ A. Salnikov'^s^ 

J. Salt^^^, B.M. SalvachuaFerrando^ D. Salvatore^^"'^^'', F Salvatorei'*^ A. Salvucci'*^ A. Salzburger^^, 

D. Sampsonidis'^^ B.H. Samset"^, H. Sandaker'^ H.G. Sander^', M.R Sanders^^ M. Sandhoff"^, 
R Sandhu'5^, R. Sandstroem"'^ S. Sandvoss'", D.RC. Sankey'^?^ g Sanny'^^ A. Sansoni"*^, 

C. Santamarina Rios*^\ C. Santoni^^ R. Santonico'^^^ '''^'', J.G. Saraiva'^"^", T. Sarangi"\ 

E. Sarkisyan-Grinbaum^ F. Sarri'22a,i22b^ q Sasaki^*', N. Sasao''^ I. Satsounkevitch^*', G. Sauvage"^, 
R Savardi57.^ A.Y. Savine^ V. Savinoyi^s, L. Sawyer^"*./, D.H. Saxon", L.R Says", C. Sbarrai^'^'i^'', 
A. Sbrizzi'^^' i^'', D.A. Scannicchio^^ J. Schaarschmidt^^, R Schacht'^^ U. Schafer**\ S. Schaetzel^^'', 
A.C. Schaffer'i^ D. Schaile'^^ R.D. Schamberger'^^, a.G. Schamovio\ V.A. Schegelsky'^i, 
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Abstract. The ATLAS detector at the Large Hadron Collider has collected several hundred 
million cosmic ray events during 2008 and 2009. These data were used to commission the Muon 
Spectrometer and to study the performance of the trigger and tracking chambers, their align- 
ment, the detector control system, the data acquisition and the analysis programs. We present 
the performance in the relevant parameters that determine the quality of the muon measurement. 
We discuss the single element efficiency, resolution and noise rates, the calibration method of 
the detector response and of the alignment system, the track reconstruction efficiency and the 
momentum measurement. The results show that the detector is close to the design performance 
and that the Muon Spectrometer is ready to detect muons produced in high energy proton-proton 
collisions. 
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1 The ATLAS Muon Spectrometer 

The ATLAS Muon Spectrometer (MS in the follow- 
ing) is designed to provide a standalone measure- 
ment of the muon momentum with an uncertainty in 
the transverse momentum varying from 3% at 100 
GeV to about 10% at 1 TeV, and to provide a trig- 
ger for muons with varying transverse momentum 
thresholds down to a few GeV. A detailed description 
of the muon spectrometer and of its expected perfor- 
mance can be found in [ 1 1 1 2 1 1 3 1 . Here only a brief 
overview is given. The muon momentum is deter- 
mined by measuring the track curvature in a toroidal 
magnetic field. The muon trajectory is always nor- 
mal to the main component of the magnetic field so 
that the transverse momentum resolution is roughly 
independent of rj over the whole acceptance. The 
magnetic field is provided by three toroids, one in 
the "barrel" (|t7| < 1.1) and one for each "end-cap" 



(1.1 < I T] I < 2.7), with a field integral between 2 and 
8 Tm. The muon curvature is measured by means 
of three precision chamber stations positioned along 
its trajectory. In order to meet the required precision 
each muon station should provide a measurement on 
the muon trajectory with an accuracy of 50 /xm. In 
Figure[T]a schematic view of the muon spectrometer 
Qis given. 

For most of the acceptance Monitored Drift Tube 
(MDT) chambers are deployed [2|. The coordinate 
in the plane perpendicular to the wires, measured by 
the MDT, is refeiTed to as the precision, or bending 
coordinate, being mainly perpendicular to the direc- 
tion of the toroidal field. In the end-cap inner region, 
for |t7| < 2.0, Cathode Strip Chambers (CSC) |2| are 
used because of their capability to cope with higher 
background rates. 

The MDT chambers are composed of two Mul- 
tiLayers (ML) made of three or four layers of tubes. 
Each tube is 30 mm in diameter and has an anode 
wire of 50 jtim diameter. The gas mixture used is 
93% Ar and 7% CO2 with a small admixture of wa- 
ter vapour, the drift velocity is not saturated and the 
total drift time is about 700 ns. The space resolution 



' The ATLAS reference system is a cartesian right- 
handed coordinate system, with the nominal collision point 
at the origin. The positive x-axis is defined as pointing 
from the collision point to the centre of the LHC ring and 
the positive y-axis points upwards while the z-axis is tan- 
gent to the beam direction at the collision point. The az- 
imuthal angle (p is measured around the beam axis, and 
the polar angle is the angle measured with respect to the 
z-axis. The pseudorapidity is defined as r7=-ln tan(}/2. 
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ATLAS 




Fig. 1. Schematic view of the muon spectrometer in the x-y (top) and z-y (bottom) projections. Inner, Middle and Outer 
chamber stations are denoted BI, BM, BO in the barrel and EI, EM, EO in the end-cap. 
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attainable with a single tube is about 80/xm, mea- 
sured in a test beam H, 131. The CSC chambers are 
multiwire proportional chambers with cathode strip 
read out. The cathode planes are equipped with or- 
thogonal strips and the precision coordinate is ob- 
tained measuring the charge induced on the strips 
making the charge interpolation between neighbour- 
ing strips. Typical resolution obtained with this read- 
out scheme is about 50 jJ-m. 

The trigger system of the MS is based on two dif- 
ferent chamber technologies: Resistive Plate Cham- 
bers (RPC) [2] instrument the barrel region while 
Thin Gap Chambers (TGC) |2| are used in the 
higher background environment of the end-cap re- 
gions. Two RPC chambers are attached to the mid- 
dle barrel chambers providing a low-pj trigger. A 
high-pT trigger is provided by the RPC modules in- 
stalled on the outer barrel chambers in combination 
with the low pr signal provided by the middle cham- 
bers. The RPCs also provide the coordinate along the 
MDT wires that is not measured by the MDT cham- 
bers. 

Similarly in the end-cap two TGC doublets and 
one triplet are installed close to the middle station 
and provide the low-pr and high-y^j- trigger signals. 
The TGCs also measure the coordinate of the muons 
in the direction parallel to the MDT wires. This co- 
ordinate is referred to as the second, or non-bending 
coordinate. For this purpose TGC chambers are also 
installed close to the MDTs in the inner layer of the 
end-cap (EI). 

Some MS naming conventions adopted in this 
paper are introduced here. The MS is divided in the 
x-y-plane (also referred to as ^ -plane) in 16 sectors: 
Sector 5 being the uppermost and Sector 13 the lower 
most. In both barrel and end-cap regions the MS is 
divided into 8 'Large' sectors (odd numbered sec- 
tors) and 8 'Small' sectors (even numbered sectors), 
determined by their coverage in (j). The muon sta- 
tions are named 'Inner', 'Middle', and 'Outer', ac- 
cording to the distance from the Interaction Point 
(IP). The three stations for the barrel are denoted BI, 
BM, and BO, and for the End-Cap EI, EM, and EO, 
respectively. Along the z axis, the MS is divided into 
two sides, called side A (positive z) and C (negative 
z). 

As a complementary source of information, two 
publications |4|, |5 1 on a detector system-test with a 
high momentum muon beam can be consulted. 



Beginning in September 2008 the ATLAS detec- 
tor was operated continuously up to November 2008 
and then for different periods starting from Spring 
2009. The first beams were circulated in the LHC 
machine in September 2008 but no beam-beam col- 
lisions were delivered. During these periods, the AT- 
LAS detector collected mainly cosmic ray data. All 
muon detector technologies were included in the run 
with the exception of CSCs for which the Read Out 
chain was still not yet commissioned and therefore 
they are not included in the results presented in this 
paper. 

The analyzed data samples and the reconstruc- 
tion software are described in Section |2] The cos- 
mic ray trigger is described in Section [3] Studies of 
data quality, calibration, and alignment are presented 
in Sections |4] |5] |6] |7] respectively, while studies on 
tracking performance are presented in sections|8]and|9] 
The results are summarized in SectionfTol 

2 Data sample and reconstruction 
software 

2.1 Data sample 

In preparation for LHC collisions, the ATLAS de- 
tector has acquired several hundred million cosmic 
ray events during several run periods in 2008 and 
2009. The analysis of a subset of data corresponding 
to about 60 M events is presented here. These runs 
allowed commissioning the ATLAS experiment, the 
trigger, the data acquisition, the various detectors 
and the reconstruction software. Most of the cosmic 
rays reach the underground detectors via the two big 
shafts. They have incident angles close to the vertical 
axis and they are mainly triggered by the RPCs. The 
selected runs, together with the status of the mag- 
netic field in the MS and the number of collected 
events for the different trigger streams, are listed in 
Table □ 

2.2 Muon reconstruction software 

The data were processed using the complete ATLAS 
software chain f6l: data decoding, data preparation 
(which includes calibration and alignment), and track 
reconstruction. Muon reconstruction has been han- 
dled by two independent packages, namely Moore flT] 
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KUIl 


Trigger 


D-ilClU 


Ui CVIS 


Period 


91060 


RPC 


Off 


17 M 


Fall 08 


91060 


TGC 


Off 


0.2 M 


Fall 08 


89106 


TGC 


Off 


0.4 M 


Fall 08 


89403 


TGC 


Off 


0.4 M 


Fall 08 


91803 


TGC 


On 


50 K 


Fall 08 


91890 


RPC 


On 


16 M 


Fall 08 


113860 


RPC 


Off 


6M 


Spring 09 


121080 


RPC 


On 


21 M 


Summ 09 



Table 1. List of analyzed data runs together with the cor- 
responding trigger stream, statistics and status of the MS 
magnetic field. All runs were collected in Fall 2008, with 
the exceptions of run 1 13860 collected in Spring 2009, and 
run 121080 in Summer 2009. 



and Muonboy [ 8 1 . The two reconstruction algorithms 
are similar in design but differ in some details. The 
general strategy is to reconstruct muon trajectories 
both at the local (individual chamber), as well as at 
the global (spectrometer), level. The trajectories re- 
constructed in individual chambers can be approxi- 
mated as straight lines over a short distance where 
bending has little effect and are therefore fit to track 
segments. Full tracks are formed by combining seg- 
ments from multiple chambers. 

Prompt muons produced in proton-proton colli- 
sions have trajectories that point back to the Interac- 
tion Point (IP). Moreover they are synchronous with 
the collision since all the detector front end electron- 
ics are synchronized with the LHC bunch crossing 
frequency of 40 MHz. In contrast, cosmic ray muons 
are "non-pointing" and are asynchronous with the 
detector clock: they have an additional 25 ns jitter 
with respect to the clock selected by the trigger In 
addition, during commissioning the different trigger 
detectors were not timed with sufficient precision, 
leading to variations in timing depending on the re- 
gion of the detector that originated the trigger. A fur- 
ther difficulty in track reconstruction was due to the 
lack of precise alignment of the muon detectors dur- 
ing this commissioning phase, as described in Sec- 
tion [T] 

The reconstruction algorithms were adapted for 
these "cosmic ray" conditions as described below. 
Both programs were modified by relaxing the stan- 
dard tracking requirements and implementing a pro- 
cedure to accommodate the cosmic ray timing condi- 
tions. The tolerance for hit association to form track 



segments and the uncertainty associated with each 
hit position were increased. Moreover, a procedure 
called global to refit (GfQ-refit) was developed in both 
reconstruction algorithms to compensate for the 25 
ns time jitter and the imprecise trigger timing. The 
aim of this procedure is to determine with better pre- 
cision the time when the cosmic ray crossed the de- 
tector by introducing a free global timing parameter 
(gfo) in the segment reconstruction. The implemen- 
tation of the Gfo-refit in the two reconstruction al- 
gorithms is briefly described below while the results 
are presented in Section|5] 



2.3 Muonboy track reconstruction 

The strategy of the Muonboy reconstruction algo- 
rithm can be summarized in four main steps: 

- identification of Regions Of Activity (ROA) in 
the muon system with the information provided 
by the RPC/TGC detectors; 

- reconstruction of local segments in each muon 
station in the identified ROA; 

- combination of segments of different muon sta- 
tions to form muon track candidates using three- 
dimensional tracking in magnetic field; 

- global track fit of the muon track candidates through 
the full system using individual hit information. 

The topology of cosmic ray tracks is accommodated 
by relaxing the Region Of Activity requirement of 
pointing in a projective geometry when associating 
hits to form segments, or matching segments to form 
tracks. Moreover, since cosmic ray events have low 
occupancy, looser quality criteria were used for the 
selection of segments and tracks. 

The Muonboy algorithm for the Gfo-refit consists 
of a scan of different gfo values in steps of 10 ns, 
doing the full segment reconstruction at each step. 
The gfo value giving the best reconstruction quality 
factor is kept and a parabolic fit is performed using 
this best value and the two closer values along the 
parabola. Then the gfo corresponding to the mini- 
mum of the quality factor parabola is chosen. In or- 
der to obtain high efficiency, the accuracy require- 
ment for the MDT single hit resolution is relaxed by 
adding in quadrature a 0.5 mm constant smearing to 
the intrinsic resolution function (described in Sec- 
tion |6]l. This smearing is increased by additional 0.5 
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mm if the Gfo-refit fails. Moreover, a less demand- 
ing track quality factor is required for tracks when 
hits are missing or are not associated to the track. 

2.4 Moore track reconstruction 

The Moore reconstruction algorithm is built out of 
several distinct stages: 

- identification of global roads throughout the en- 
tire spectrometer using all muon detectors (MDT, 
CSC, RPC and TGC); 

- reconstruction of local segments in each muon 
station seeded by the identified global roads; 

- combination of segments of different muon sta- 
tions to form muon track candidates; 

- global, three-dimensional, tracking and final 
track fit. 

Several modifications to the standard pattern 
recognition were made to optimize the reconstruc- 
tion of cosmic ray tracks. In the global road finding 
step, a straight line Hough transform was used to al- 
low for non-pointing tracks. The cuts on distance and 
direction between the road and the segment were re- 
laxed. In the segment finding no cuts were applied 
on the number of missing hits {i.e. drift tubes that 
are expected to be crossed but have no hits). 

The Gfo-refit consists in varying simultaneously 
the global time offset (gfo) for each segment recon- 
structed in a chamber Then all measured times of 
hits associated to the segment are translated into drift 
radii after subtraction of the gfo. The gfo value that 
minimizes the sum in quadrature of the weighted 
residuals (corresponding to the segment reconstruc- 
tion x^) is selected. 

In this fit the MDT uncertainties are set to twice 
the test-beam drift tube resolution. If the segment fit 
is not successful, a straight line fit is performed as- 
suming a constant 1 mm error. Hits are removed if 
their distance from the segment is greater than la. In 
the track fit the MDT errors are enlarged to 2 mm to 
account for uncertainties in the alignment of cham- 
ber stations. 

3 Trigger configuration during data 
talking 

A more detailed description of the trigger system can 
be found in li2J,|i9J. Here only specific issues related 



to the 2008-2009 cosmic ray data taking are intro- 
duced. The muon level- 1 trigger is issued by the RPC 
in the barrel and by the TGC in the end-caps. During 
cosmic ray data taking most of the statistics were col- 
lected using this trigger. Special trigger configura- 
tions were adopted with different geometries (eg non 
pointing to the IP) and different timing (e.g. delay- 
ing the triggers issued by the upper sectors in order 
to trigger only in the lower sectors to mimic particles 
coming from the IP) when commissioning the muon 
trigger system itself or when selecting cosmic rays 
for commissioning the other ATLAS sub-systems. 

In beam-collision configuration, the level- 1 
muon trigger selects pointing tracks with six differ- 
ent thresholds in transverse momentum and sends 
information to the Central-Trigger-Processor (CTP). 
The six thresholds, three low-pj- and three high- 
Pt, do not distinguish between different detector re- 
gions, barrel or end-cap. For cosmic rays, to help 
commissioning separately the two regions, it was 
chosen to assign three thresholds to the barrel and 
three to the end-cap. 

3.1 Barrel level-1 trigger 

The barrel trigger detectors are arranged in three sta- 
tions each having a doublet of RPC layers at increas- 
ing distances from the IP. In each sector the first two 
stations are mechanically coupled to the BM MDT 
while the third is coupled with the BO MDT as shown 
in Figure [T] 

The trigger algorithm is steered by signals on 
the middle layers, named Pivot plane. When a hit is 
found on this plane, the I0W-/97- trigger logic searches 
for hits in the inner layers, named Confirm plane, and 
requires a coincidence in time of three hits over the 
four layers in a pre-calculated cone. The width of this 
cone defines the pj threshold. If hits are also found 
in a pre-calculated cone of the outermost plane in 
coincidence with a I0W-/97- trigger, a high-pr trig- 
ger is issued. Also in this case the pj threshold is 
defined by the width of the cone. In addition to the 
Pt requirement, the trigger logic also demands the 
track to be pointing towards the IP both in and 77 . 
In the cosmic ray runs only three of the six thresh- 
olds were used in the barrel and were defined as 
MU0±OW, MUOJiIGH and MU6. The two thresh- 
olds MUOJLOW/HIGH did not select a physical pr 
range; in fact, the MUOJ^OW was triggered only 
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by the time coincidence of 3 out of 4 hits with- 
out any pointing constraint and the MUOJ-IIGH was 
triggered by the coincidence of a MUOJ^OW with 
at least a hit in the corresponding outer plane. The 
threshold MU6 required not only a time coincidence 
but also an IP-pointing constraint in the -projection 
only. To emulate the timing expected for beam col- 
lisions, and to enhance the illumination of the In- 
ner Detector (ID), the cosmic ray trigger was issued 
mainly by the bottom sectors. This was achieved by 
delaying the top sector trigger by 5 BC (125ns) pre- 
venting it from arriving first at the Central Trigger 
Processor (CTP) and thus forming the trigger. 

In the fall 2008 data taking period, the timing 
of the low- pt trigger and the data read-out laten- 
cies were still under commissioning. This had a large 
impact on the detector coverage. The situation has 
largely improved for the runs taken in 2009 both in 
terms of detector coverage and in trigger timing as 
shown in Sections [4. 3l and l672l 

3.2 End-cap level-1 trigger 

The level-1 TGC trigger system provided three 
thresholds, named MU0.TGC.HALO, MUOJTGC 
and MU6-TGC. The trigger was issued by the co- 
incidence between several TGC layers. The logic 
was based both on timing (BC identification) and 
geometry (pointing track). The main difference be- 
tween the three trigger thresholds is related to the 
required number of layers and to the degree of point- 
ing to the IP MUO.TGCJIALO required a 3 out of 4 
layer coincidence in the two outermost TGC stations, 
the so-called Doublet chambers, in both rj (bending) 
and (p (non-bending) projections and a pointing re- 
quirement within 20°. MUO.TGC and MU6.TGCie- 
quired in addition a 2 out of 3 layer coincidence in 
the TGC stations closer to the IP, the so-called Triplet 
chambers, in the rj projection only. The pointing re- 
quirement of MUO-TGC was of ±10° degrees while 
for MU6.TGC was of ±5°. 

The trigger was timed for high-momentum 
muons coming from the IP. All the delays due to 
different time-of-flight and cable lengths were prop- 
erly set and cross-checked using a test pulse system 
achieving a relative timing within 4 ns. For most of 
the cosmic run period, only the level-1 trigger gen- 
erated from the TGC bottom sectors was used. This 
was chosen to ensure good timing of the trigger with 



the read-out of the ID, since cosmic muons triggered 
by the TGC bottom sectors and crossing the ID have 
a time-of-flight similar to muons produced in colli- 
sions. 

4 Data quality assessment 

4.1 Introduction 

The data quality assessment consists of several soft- 
ware algorithms working at different levels of the 
data taking. The Detector Conti'ol System (DCS) HIl 
is the first source of information available during the 
operation of the detector Here information on the 
hardware status of the different sub-detectors and on 
the settings of Low Voltage (LV) and High Voltage 
(HV) power supplies and on the gas system is avail- 
able. The DCS also receives information from the 
Data Acquisition (DAQ) as soon as problems 
during the read-out of a chamber appear 

The next stage in the chain of data quality as- 
sessment is the on-line monitoring. It receives in- 
put from the data acquisition system running in a 
spectator mode. Once the data are decoded, monitor- 
ing histograms are filled showing quantities related 
to the detector operation. Part of the muon data se- 
lected by the level- 1 trigger Region Of Interest (ROI) 
are transferred by the level-2 trigger processors to 
three dedicated computing farms (referred to as cal- 
ibration centers) to monitor and determine the cali- 
bration parameters of the MS chambers. The larger 
event samples available at the calibration centers al- 
low the analysis of single drift tube responses. The 
goal of the analysis at the calibration centers is to 
provide drift tube and trigger chambers calibration 
constants and to give general feedback on the de- 
tector operation within 24 hours, which is the time 
needed, at high luminosity, to collect enough statis- 
tics to calculate new calibration constants. 

On a longer time scale, using the full reconstructed 
ATLAS event information, the off-line data monitor- 
ing provides the final information on the data quality. 
At each step a flag summarizing the data quality at 
that level is stored in a database. 

4.2 IVIDT chambers 

In the fan 2008 period (e.g. Run 91060) only five 
out of 1110 MDT chambers were not included in the 
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data taking. Of these five chambers, two were not yet 
connected to services and three had problems with 
the gas system. Due to the cosmic ray illumination 
and the trigger coverage not all chambers had suf- 
ficient event samples to determine the performance 
of single drift tubes. The studies reported here were 
done at different levels of detail, from chamber infor- 
mation down to single drift tube information when 
the event samples were sufficient. The data survey 
searched for problems of individual read-out chan- 
nels as well as of clusters corresponding to hard- 
ware related groups of tubes. A screen shot of one 
of the online monitoring applications used for the 
MDT chambers is shown in Figure |2] Here the av- 
erage number of hits per tube for each MDT is rep- 
resented in a T]-^ plot where the higher cosmic il- 
lumination on the top and bottom sectors (3-7, 11- 
15) compared to the vertical sectors (16- 2 and 8-10) 
is clearly seen as well as the larger illumination on 
the A side of the detector where the larger shaft is 
present. The five chambers not included in the data 
acquisition are marked as dark gray boxes. Two more 
chambers are visible with very low statistics due to 
problems with the HV supplies. For 32 MDT cham- 
bers one of the two multi-layers was disconnected 
from HV. 

A detailed list of hardware problems found in run 
91060 is reported in Table|2] The cosmic ray flux was 
not sufficient for a detailed analysis of single drift 
tubes for 15 MDT chambers (^3K channels). Thus 
we were able to analyze individually 336K, out of 
the working 339K, drift tubes. To summarize, about 
5K channels, out of 336K, have shown some prob- 
lems in run 91060, corresponding to 1.5%. Most of 
these channels have been recovered during the 2008- 
2009 shutdown period. Only a very small fraction of 
problems, at the level of a few per mill, could not 
be solved, such as permanently disconnected tubes 
(broken wires) or chambers with very difficult ac- 
cess. 

In addition to monitoring in the DAQ framework 
(on-line monitoring), the data are also processed with 
the offline reconstruction program which produces 
monitoring histograms. This ensures that the recon- 
struction works properly and that the correct condi- 
tions data (calibration and alignment constants) are 
used in the first processing of the data. The off-line 
monitoring gathers and presents information on sev- 
eral variables for single drift tubes, e.g. drift time 



and collected charge distributions, hit occupancy and 
noise rate. These variables are obtained for individ- 
ual MDTs or grouped for regions, such as 7] or 
sectors, barrel or end-cap, side A or C. Variables re- 
lated to segments or tracks are also monitored. 



4.3 Barrel trigger chambers: RPC 

Commissioning of the RPC detectors progressed 
continuously and substantial improvements were 
made during the 2008-2009 shut-down. As an ex- 
ample Figure |3] shows a two-dimensional distribu- 
tion of RPC strips requiring a 3 out of 4 majority co- 
incidence for the low- pt trigger demonstrating that 
the trigger coverage in Spring 2009 was at the 95 % 
level. 

Studies of the trigger performance were made us- 
ing the data of run 91060 after implementation of the 
trigger roads |10|. For the low-pj trigger the four 
RPC layers in the Middle station are involved, both 
in Tj and (j) projections. Tracks were accepted by the 
trigger if any strip of the pivot plane was in coin- 
cidence with a group of strips of the confirm plane 
aligned with the IP, realizing a majority combination 
of 3 out of 4 RPC layers. Figure|4]left shows the spa- 
tial correlation between (p strips in the pivot planes 
and (j) strips in the confirm planes. The correlation 
line is slightly rotated with respect to the diagonal 
due to the different distance of the confirm and pivot 
planes with respect to the IP. 

A random trigger was used to measure the count- 
ing rate for each read-out strip. This is a measure- 
ment of the RPC system noise rate. About 3 lOK strips 
were analyzed over a total of 350K working strips. 
Figure |4]right shows the distribution of single chan- 
nel noise rate, normalized to an area of 1 cm^. For 
each strip, the noise rate is calculated as the number 
of hits divided by the number of random triggers and 
the width of the read-out gate of 200 ns, and is nor- 
malized to the area of the strip (typically 550 cm^ 
for a BM eta strips and 900 cm^ for a BO eta strips). 
Only a few hundred strips showed a counting rate 
above 10 Hz/cm^ which is the background rate ex- 
pected when the LHC will run at high-luminosity. 
The average noise rate of the RPC was stable during 
the different running periods. 

The fraction of dead channels, considering only 
the part of the detector included in the read-out in the 
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Fig. 2. Screen shot of a monitoring application displaying the MDT hit occupancy for all chambers. Each chamber is 
represented by a small box. The color of the box is related to the average number of raw hits per tube. The boxes are 
arranged in an rj-tj) grid: a column represents an Tj slice, perpendicular to the beam axis; a row represents one of the 
sixteen (j) sectors. Within each sector chambers of the Inner, Middle, Outer ring are displayed separately. 



Number of channels analyzed with sufficient event samples 


336144 


Fraction 


Channels not included in the read-out 
Channels with read-out or initialization problems 
Channels with HV or gas problems 
Permanently dead channels (broken wires) 


936 
744 
2942 
323 


0.28% 
0.22% 
0.88% 
0.10% 


Total problematic channels 


4945 


1.47% 



Table 2. List of MDT channels with problems in run 91060. 



Fall 2008 runs, was 1.5%, mainly due to problems in 
the front-end electronics. 



4.4 End-cap trigger chambers: TGC 

In the end-caps the muon trigger is provided by the 
TGC chambers installed in three layers that surround 
the MDT Middle chambers. All together they form 



the so-called Big Wheels (BW), one in each end-cap. 
In addition, TGC chambers are also installed close to 
the EI chambers in the Small Wheels (S W), but these 
are only used to measure the muon (p coordinate. In 
Fall 2008 all the BW TGC sectors were read-out. 
Given the installation schedule for the ATLAS de- 
tectors, the Inner TGC station were the last chambers 
installed and they were not fully operational during 
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Fig. 3. RPC \ow-pj trigger coverage in rj-tj) for Run 
1 13860 (spring 2009). Eacli rj and (p strip producing a low- 
pr trigger corresponds to an entry in tlie plot. The cover- 
age in spring 2009 was about 95%. 



2008 runs. For this reason they are not discussed in 
the following. 

Two types of trigger configuration were adopted 
in Fall 2008. One was optimized to study the end- 
cap muon detectors with cosmic rays. In this config- 
uration all TGC BW sectors were used in the trig- 
ger. The other setting was optimized to provide the 
trigger for the ID tracking detectors and was used 
for timing the ID. In order to mimic muons com- 
ing from the IP, only the five bottom sectors were 
used to trigger. The typical detector coverage in these 
two trigger configurations is shown in Figures |5] by 
plotting the coincidence positions in the x-y plane 
for wire and strip hits for run 91060 (left) and run 
91803 (right). Only about 0.8% of chambers were 
not operational due to HV or gas problems. Since 
for the trigger a majority logic is required these in- 
active chambers do not produce any dead regions in 
the trigger acceptance. 

The HV and front-end threshold setting, the gate 
widths for wires and strips, and the trigger sectors 
are listed in Table[3]for these two runs. 

For each trigger issued by the CTP, the TGC Read 
Out Driver (ROD) sends to the DAQ system the data 
corresponding to three Bunch Crossings (Previous, 
CuiTent and Next BC) contained in two separate buffers. 
Of the two buffers, one is located in the front-end 
board where the wires and strips providing the low- 
Pt coincidence are separately recorded. In the sec- 
ond buffer, located in the Sector Logic Board in the 
service counting room, the coincidence of the wire 
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Fig. 4. A): RPC spatial correlation between the pivot strip 
number and the confirm strip number in the (j) projection 
for a programmed trigger road. 128 strips correspond to a 
RPC plane 3.8 m long. B): distribution of strip noise rates 
per unit area measured with a random trigger for 31 OK 
RPC strips. The larger noise present on some strips is prob- 
ably due to local weaknesses of grounding connections. 



and strip signals is done. Each buffer has a programmable 
identifier that has to be adjusted in order to read out 
the coiTect (Current) BC data. Figure |6] shows the 
readout timing for the front-end and the sector logic 
buffers for level- 1 triggers issued by the TGC. About 
98.6% of data in the front-end buffer, and 99.8% of 
data in the sector logic buffer are read out with the 
correct timing. The small population in the previous 
or next BC is due to cosmic ray showers. 
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Run 


Trigger sector 


HV 


Threshold 


Gate widths for wire / strip 


91060 


8 to 12 


2800 V 


100 mV 


35 /45 ns 


91803 


1 to 12 


2650 V 


80 mV 


35 / 45 ns 



Table 3. TGC sectors participating in the trigger, high voltage setting, threshold and gate width. 
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Fig. 5. Map of coincidences of wire and strip hits in 
the x-y plane. A): the five bottom sectors (sectors 8-12, 
195° <if < 345°) used for timing the ID tracking detec- 
tors in run 91060. B): with all sectors participating in the 
trigger during run 91803. 
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Fig. 6. A): TGC front-end and B): sector logic buffers 
for BC identification. Three BC crossing, previous, current 
and next are readout. 



5 MDT chamber calibration 
5.1 Calibration method 

The MDTs require a calibration procedure lfT2l to 
precisely convert the measured drift time into a drift 
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distances from the anode wire (drift radius) that is 
subsequently used in pattern recognition and track 
fitting. The calibration of the MDT chambers is per- 
formed in three steps. In the first step the time offset 
with respect to the trigger signal, fo, for each tube 
or group of tubes is determined; in the second step 
the drift-time to space relation, r{t) function, is com- 
puted; in the third step the spatial resolution is deter- 
mined. 

The calibration constants are loaded in the Con- 
ditions Data Base (known as 'COOL') [13] and then 
retrieved, according to an Interval Of Validity (lOV) 
mechanism, to be used in the offline reconstruction. 
The lOV determines for which group of runs the cal- 
ibration constants are valid. The gas mixture com- 
position varied during the data taking period since 
the water injection part of the gas system was un- 
der commissioning resulting in a not constant ad- 
mixture of water vapour, as can be seen in Figure |9] 
Nonetheless the calibration procedure based on the 
lOV mechanism was able to provide good calibra- 
tion constants for all the running period. 

The fo offset depends on many fixed delays like 
cable lengths, front-end electronics response, Level- 
1 trigger latency, time of flight from the IP and has 
to be determined for each drift tube. The offset is 
obtained by fitting a Fermi function to the leading 
edge of the drift time distribution as shown in Fig- 
ure |7] left. The precision expected in LHC collision 
data is better than 1 ns with a dataset of about lOK 
muons crossing the drift tube. This uncertainty does 
not significantly degrade the position resolution of 
the MDT tubes which corresponds to a time span of 
about 5 ns. In Figure|2]right also the typical spectrum 
of ADC for all tubes in a chamber is reported. Charge 
information in each tube is obtained using a Wilkin- 
son ADC 1 14 |. As the MDT chambers are operated at 
different temperatures depending on their positions 
in the MS, the r{t) functions differ depending on 
location and are determined separately. In addition, 
variations of the toroidal magnetic field along the 
drift tubes produce different Lorentz angles, thus dif- 
ferent r{t) functions. An initial rough estimate of the 
r(f ) function is obtained with an accuracy of 0.5 mm 
by integrating the drift-time distribution. This is cor- 
rect under the approximation of a uniform (i«/iir dis- 
tribution, where n is the number of hits at a drift ra- 
dius r 




ADC counts 

Fig. 7. A): typical drift time spectrum in cosmic ray events 
for an MDT chamber. The position of the inflection point 
of the leading edge of the spectrum, to, is determined by 
fitting a Fermi function (shown in red) to the beginning of 
the spectrum. B): typical spectrum of ADC for all tubes 
in a chamber. Hits below 50 ADC counts are identified as 
electronic noise. 
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Nhits is the total number of hits in the time spectrum 
and r,„ux is the maximum drift radius (14.4 mm). In 
cosmic rays this approximation is only good at the 
level of a few hundred |j.m mainly because of the 
production of 5 -ray electrons along the track. An 
r(f ) relation with a higher accuracy, of about 20 |J.m, 
is obtained from this initial estimate by applying cor- 
rections, 5r{t), which minimize the residuals of track 
segment fits with an iterative procedure. This min- 
imization procedure, called auto -calibration, takes 
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into account the dependence of the parameters of the 
fitted segments on the applied corrections 5r{t) and 
is mainly based on the geometrical constraints from 
the precise knowledge of the wire positions. Figure[8] 
shows a typical residual distribution of a chamber, as 
a function of the distance of the track segment from 
the anode wire, after the auto-calibration. 

In cosmic ray events additional sources of time 
jitter, beyond the intrinsic resolution, spoil the MDT 
measurement. The first cause of time jitter is due to 
cosmic ray muons crossing the tubes with an arbi- 
trary phase with respect to the front-end electronics 
clock |15|. This implies a time jitter corresponding 
to a 25 ns uniform distribution. The second cause 
is related to the spread of the trigger time for trig- 
gers generated in different parts of the detector (up to 
about 100 ns due to the initial stage of the trigger tim- 
ing). Two different methods have been alternatively 
used to reduce the impact of these effects: the RPC- 
time correction and the MDT Gt^-refit. The achieved 
performance with both methods are discussed in Sec- 
tion|6] In the following a brief description of the for- 
mer method is given. 

The RPC-time correction uses the trigger time 
measured by the RPC chambers on an event by event 
basis. This time correction was applied only to the 
MDT chambers of the BM stations since these cham- 
bers are close to the two RPC stations used to is- 
sue the trigger and so no corrections due to time 
of flight and, more importantly, no corrections due 
to the spread in timing of the trigger signals issued 
by different parts of the detector are needed. This 
method cannot be applied to the end-cap region since 
the TGC do not provide a measurement of the trigger 
time but rather they select the appropriate BC. 

With this correction the time jitter due to the two 
effects mentioned above is reduced from ^100 ns 
to few ns (see Section |6]l. The distribution of the 
residuals obtained after calibration using the RPC- 
time correction method is presented in Figure[8]left. 
The precision of the auto-calibration is better than 
~20 i^m using this correction. 

The Gfo-refit has also been used to improve the 
single tube resolution, as discussed in Section|6] Also 
the precision of the auto-calibration is much improved 
with respect to the uncorrected situation. As shown 
in Figure |8] a precision of ~'50 \iva is obtained for 
the residuals of the segment fit after auto-calibration 
with small residual systematics on the auto-calibration. 




Fig. 8. A): Residuals as a function of the track segment 
distance from the wire after the r{t) auto-calibration and 
RPC-time corrections. The points correspond to the mean 
value of the distribution of residuals and the error bars to 
its RMS value. B): Residuals as a function of the track seg- 
ment distance from the wire after the r{t) auto-calibration 
using the Gfo-refit method. The points correspond to the 
mean value of the distribution of residuals and the error 
bars to its RMS value. Residual systematics at the level of 
50 |J.m are present using this correction. 
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5.2 End-cap chambers calibration with 
monitoring chamber 



6 Detector performance: efficiency 
and resolution 



For the end-cap MDTsystem, due to the limited num- 
ber of cosmic ray events, a different method to deter- 
mine the r(f ) relation was used. A small MDT cham- 
ber installed on the surface of the ATLAS under- 
ground hall was set up 1 16| to monitor continuously 
the MDT gas composition. One multi-layer is con- 
nected to the supply line of the gas recycling system 
while the other is connected to the return line. This 
chamber benefits from a very large cosmic ray rate 
and can therefore determine the r(f) function with 
high precision in short time intervals. Cosmic ray 
muons are triggered by scintillator counters mounted 
on the monitoring chamber. The trigger time is mea- 
sured and subtracted event-by-event from the tube 
drift times, in this way the jitter related to the asyn- 
chronous front-end clock is automatically removed. 

Data from the monitoring chamber are used to 
derive a r(f ) function every 6 hours to monitor the 
gas drift properties. Figure |9] shows the variation of 
the maximum drift time (the drift time of muons cross- 
ing the drift tube close to its edge) over the period 
September-October 2008. Two r{t) functions were 
used to cover the Fall 2008 run period, for each pe- 
riod the r{t) function for each chamber of the MS 
was corrected to account for the temperature differ- 
ence using the data measured by the sensors mounted 
on any particular chamber. The temperature correc- 
tions to the r(f ) function were derived from the Garfield 
MagBoltz simulation program fnl, fH,!!!). The 
output of the simulation was validated by several mea- 
surements with a muon beam |5|. In the end-cap re- 
gion, the temperature varies by about 4°C from top 
to bottom of the MS, resulting in a variation of the 
maximum drift time of about 10 ns, nn the other hand 
the temperature of the cavern was remarkably stable 
in time. 

Figure[TO]shows the distribution of the mean and 
RMS value of the residuals from the fit to track seg- 
ments in all end-cap chambers (run 91060). A Gaus- 
sian fit is superimposed. The r(f) function derived 
from the gas monitor chamber with temperature cor- 
rections provides an acceptable calibration for all the 
MDT chambers of the end-cap: the average standard 
deviation of the residuals is about 100 |J.m. 



6.1 IVIDT 

6.1.1 MDT drift time distribution 

The behaviour of the drift time distributions of in- 
dividual tubes is an important quality criterion for 
the MDT performance. The minimum and maximum 
drift times, to and t,nax, respectively, correspond to 
particles passing very close to the wire and close to 
the tube walls, and their stability indicates the stabil- 
ity of the calibration. The number of hits recorded 
in a small time window before the rising edge of 
the drift time distribution fo can be used to estimate 
the rate of noise due to hits not correlated with the 
trigger. A precise knowledge of to for each tube is 
essential for high quality segment and track recon- 
struction. As explained in Section|5] for cosmic rays 
some additional time jitter is present and must be ac- 
counted for 

In order to improve the quality of track recon- 
struction the Gfo-refit time correction has been used. 
The performance of the Gfo-refit algorithm has been 
investigated in the past, both using simulated data 
and using data taken with a BIL (Barrel Inner Large) 
chamber in a cosmic ray test stand under controlled 
trigger conditions ||20l . The achieved Gfo resolution 
ranged between 2 and 4 ns depending on the cham- 
ber geometry (8 layer chambers have better resolu- 
tion than 6 layer chambers) and hit topology. In par- 
ticular the Gfo-refit algorithm cannot work if all the 
hits are on the same side of the wires, typically for 
tracks at 30° with respect to the chamber plane. The 
selection of good quality segments requires a mini- 
mum of five MDT hits and segments with all hits on 
the same side of the wires are removed. 

In addition to the Gfo-refit also the RPC-time 
correction method was used for the MDT chambers 
in the middle barrel station (BM) which are located 
closely to the RFC trigger chambers. The time mea- 
sured by these RFC can be used to correct for a global 
time offset. An example of the effectiveness of the 
method is given in Figure |7] where the drift time dis- 
tribution for a BML chamber is shown after RPC- 
time corrections. The steepness of the rising edge, 
measured as one of the parameters of the Fermi dis- 
tribution, is improved, passing from 22 ns without 
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Fig. 9. Maximum drift time measured by tlie gas monitor cliamber versus time during September-October 2008. The 
red points refers to the return line and the blue points to the supply line (green and light blue points are the time average 
of the supply and return line measurements). The large variation seen between middle of September and 10''' of October 
is due to the change of the quantity of water vapour added to the standard mixture. 



correction, to 3 ns with RPC time corrections, a value 
in agreement with results from muon beam tests |5 1. 
The precision of the RPC-time correction is about 
2 ns as explained in Section 16.21 This also includes 
the contribution of the signal propagation time in the 
RPC strips. 



The distribution of the difference between the fit- 
ted GtQ and the RPC timing correction per segment 
is shown in Figure[TT]for a BML chamber. The stan- 
dard deviation of about 4 ns is consistent with an un- 
certainty of 2 ns from the RPC-time correction added 
to an uncertainty of 3 ns introduced by the Gfo-refit 
method. Tails up to 30 ns are present in the distribu- 
tion due to bad hit topologies and background hits. 



6.1 .2 Drift tube spatial resolution 

The MDT single tube resolution, as a function of 
drift distance, was studied using different time cor- 
rections. The extraction of the resolution function is 
based on an iterative method. At the first iteration 
an approximate input resolution function is assumed. 
Only segments with a minimum of six hits are con- 
sidered. These segments are fitted again after remov- 
ing one hit at the time. Subsequently, the width of 
the distribution of the residuals for the excluded hit, 
j, is computed as a function of the drift distance from 
the wire, afitj{r). The errors of the straight line fit 
(depending on the assumed tube resolution) are then 
propagated to the excluded hit. The resolution a, (r) 
is then computed by quadratic ally subtracting from 
the standard deviation of the residuals the fit extrap- 
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Fig. 10. A) Distribution of the RMS and B) of the mean 
values of the residuals from the fit to track segments in 373 
end-cap chambers using the r{t) function derived from the 
gas monitoring chamber. The black lines represent Gaus- 
sian fits to the distributions. 
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The procedure is iterated using the new resolution 
function until the input and output resolutions agree 
within statistical uncertainties; a small number of it- 
erations (two to four) is usually needed. 

In Figure [12] the tube resolution obtained for a 
BML chamber is shown as the green band. The width 
of the band accounts for the systematic uncertainty 
of the method. Also shown (solid line) is the resolu- 
tion function obtained for an MDT chamber at a high 
energy muon test beam |5 1 with well controlled trig- 
ger timing. This can be considered as reference for 
the single-tube resolution. The resolution function 
measured with cosmic rays is consistent with a time 
degradation of the reference resolution of about 3 ns. 
This is in reasonable agreement with the 2 ns time 
resolution quoted for the RPC-time correction in ad- 
dition to a small contribution from multiple scatter- 
ing and individual tube differences in fo. 
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Fig. 11. Difference between the g/o obtained with the Gfo- 
refit method and with RPC-time correction. The width of 
the distribution is a convolution of the uncertainties of the 
RPC-time correction and the G/Q-refit method. 



The single hit spatial resolution was determined 
also by applying the Gfo-refit method to track seg- 
ments reconstructed in the same chamber The pro- 
cedure was similar to that presented above with the 
convergence of the method driven by the estimate of 
the residual pulls. The resolution function is shown 
as the blue hatched band in Figure[T2l The measured 
resolution is consistent with the test beam measured 
resolution provided that an additional time uncer- 
tainty of about 2-3 ns is taken into account. 



6.1.3 Drift tube noise 

The level of noise can be measured in each drift tube 
by looking at the drift time distribution in a given in- 
terval before Jq where only hits uncorrected with the 
trigger are present. The noise rate is obtained by di- 
viding the number of hits normalized to the number 
of triggers by the chosen time interval. The charge of 
drift tube signals, at nominal running conditions, is 
well above the ADC pedestal corresponding to about 
50 counts, see Figure [T] In the reconstruction algo- 
rithms only hits with charge above this value are con- 
sidered. The distribution of noise rate with and with- 
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Fig. 12. Drift tube resolution as a function of the radius. The green shadowed (RPC correction method) and the blue 
hatched (G/o-refit method) bands represent the resolution function measured with cosmic rays with the two different 
methods described in the text. The solid line represents the resolution measured with a high momentum muon beam [5]. 



out the ADC charge cut is shown in Figure [13] for 
all MDT drift tubes. The average noise rate is only 
60 Hz without the ADC cut and 13 Hz with ADC 
cut, the former figure corresponds to an average tube 
occupancy of less than 10^^. 

6.1 .4 Drift tube efficiency 

The single tube efficiency was studied by recon- 
structing segments in a chamber using all tubes ex- 
cept the one under observation i.e. excluding one 
MDT layer at the time in segment reconstruction. 
Two different types of inefficiencies can be defined: 
i) absence of a hit in the tube; ii) a hit is present but 
is not associated to the segment because its residual 
is larger than the association cut. The inefficiency of 
type i), referred to as hardware inefficiency, is very 



small, mostly occurring at large drift distances, near 
the tube wall, where the short track length results in 
fewer primary electrons or due to the track passing 
through the dead material between adjacent tubes. 
The inefficiency of type ii), refeiTed to as tracking 
inefficiency, is dominated by 5-electrons, produced 
by the muon itself, which can mask the muon hit if 
the 5-electron has a smaller drift time than the muon. 
Tube noise can be an additional source of this type 
of inefficiency. 

Figure [14] shows the distribution of the signed 
residuals for hits in the tube of one baiTel chamber 
as a function of the distance of the segment from the 
wire. A large population at small values of the resid- 
ual, compatible with the spatial resolution, is visible. 
Large positive residuals are associated with early hits 
mainly due to 5-electrons. If a hit is not found in the 
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Fig. 13. Distribution of the drift tube noise rate witli (shadowed histogram, bottom statistical box) and without (empty 
histogram, top statistical box) the ADC cut described in the text. In the right plot the logarithmic scale allows observation 
of the very few noisy tubes. The tail of the distribution is due to very few noisy tubes that are suffering from pick up of 
high frequencies through the HV cables or interferences due to the digital clock present in the front end electronics. 



tube traversed by the muon (thus a residual cannot 
be computed) a value of 15.5 mm is assigned, larger 
than the tube radius of 15 mm. The population of 
missing hits is visible at the top of Figure [14] and it 
peaks close to the tube wall. 

The tracking efficiency is defined as the fraction 
of hits with a distance from the segment smaller than 
n times its error, this error being a convolution of 
the tube resolution and the track extrapolation un- 
certainty. Figure [15] shows the hardware efficiency 
and the tracking efficiency as a function of the drift 
radius for n = 3, 5, and 10. The tracking efficiency 
decreases with increasing radius, mainly due to the 
contribution of 5-electrons. The average tube hard- 
ware efficiency is 99.8%; the tracking efficiency is 
97.2%, 96.3% and 94.6% for n equals to 10, 5 and 3 
respectively. 

Figure [16] shows the average value of the track- 
ing efficiency for each tube of a BML chamber for n 
= 5. The average value is about 96%. An efficiency 
consistent with zero was obtained for two tubes as 
can be seen in the expanded view on the right plot. 
These were recognized as tubes with disconnected 
wires and were not considered in the average value. 

The results of a study on all the barrel cham- 
bers with enough cosmic ray illumination to allow 
the determination of the single tube efficiency is pre- 



sented in Figure [TT] The distribution of the tracking 
efficiency for a 5(7 hit association cut is shown for 
about 8 IK drift tubes. In addition to about 0.2% of 
dead channels, less than 1% of tubes have tracking 
efficiency below 90%, mainly due to calibration con- 
stants determined with insufficient precision. 



6.2 RPC 

In addition to providing the baiTel muon trigger, the 
RPC system is also used to identify the BC of the 
interaction that produced the muon. This requires a 
time resolution much better than the bunch crossing 
period of 25 ns. For this, the time of the strips that 
form the trigger coincidence is encoded in the front- 
end with a 3 -bit interpolator providing an accuracy 
of 3.125 ns ifTOl . The distribution of the time differ- 
ence between the two layers of a pivot plane in the 
projection was used to determine the RPC time res- 
olution. With this method there is no need to correct 
for the muon time of flight and the signal propaga- 
tion along the read-out strips. The RMS width of the 
distribution shown in Figure[T8]is 2.5 ns. From this a 
time resolution of 1.8 ns is derived for the two RPC 
layers forming the coincidence. For this measure- 
ment only strips associated to a reconstructed muon 
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Fig. 14. Distribution of tlie hit residuals for tubes excluded in the segment fit, as a function of the distance of the track 
from the wire. Small residuals are associated with efficient hits. The triangular region is populated by early hits produced 
by 5-electrons. Missing hits, as explained in the text, are assigned a residual value of 15.5 mm. 



track and belonging to events with one and only one 
RPC trigger were considered. 

Two other important RPC quantities related to 
the detector performance are the efficiency and the 
spatial resolution. In order to determine the RPC effi- 
ciency two main issues have to be taken into account. 
The first one is due to the fact that the RPCs are ac- 
tually providing the muon trigger thus resulting in a 
trigger bias on the efficiency calculation. The second 
one is caused by the fact that the RPC hits are also 
used in the track reconstruction; in particular, they 
measure the coordinate in the non-bending, (j) pro- 
jection. The second effect has negligible contribution 
if the efficiency is measured for the 77 strips, since 
in this projection the track reconstruction is driven 
by the MDT. For the efficiency measurement, MDT 
tracks were extrapolated to the RPC plane and the 
layer was counted as efficient if at least one rj hit 
was found with a distance of less than 7 cm from the 
extrapolation. The effect of the trigger bias has been 
removed from the efficiency measurement of an RPC 



plane by selecting all the events in which the other 
three planes (in the case of a Middle Station) were 
producing hits, since the trigger requirement is a 3 
over 4 planes majority. The distribution of the effi- 
ciency, averaged over each layer, for the RPC cham- 
bers in the Middle stations is shown in Figure [19] 
the distribution is peaked at an efficiency of 91.3%. 
To check the remaining impact of the trigger bias on 
the efficiency measurement, the same analysis was 
repeated with a sample of cosmic rays selected with 
a calorimeter trigger (Level- ICalo trigger) indepen- 
dently of the RPC trigger response. The result for 
the efficiency is superimposed in Figure [191 ^ good 
agreement between the two distributions is observed. 

The spatial resolution is related to the clusters 
size, that is the number of strips associated to a muon 
track. A muon crossing the detector near the center 
of a readout strip, will in general produce a cluster 
of size one, while clusters of size two are only ob- 
served when muons hit a narrow region at the bound- 
ary between two strips. The actual sizes of the re- 
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Fig. 15. Tube efficiency as a function of the drift distance 
averaged over all tubes of a BML chamber. Shown are 
the hardware efficiency and the tracking efficiency for hit 
residuals smaller than 3,5, and 10 times the standard devi- 
ation of the distribution. 



gions corresponding to clusters of size one and two 
depends on the detector operating parameters, but it 
is in general true that the latter is smaller than the 
former This implies that the spatial resolution must 
be smaller when measured on a subset of data with 
only clusters of size two. The spatial resolutions of Tj 
strips was determined selecting muon tracks recon- 
structed in the MDT as explained above. For each 
RPC read out plane, the distribution of the distance 
from the extrapolated track was obtained separately 
for clusters of size one and two and then was fitted 
with a Gaussian. The RMS widths of the fit were di- 
vided by the strip pitch (ranging from 27 to 32 mm 
depending on the chamber type) to allow for compar- 
ison between different RPC and are shown in Fig- 
ure |20] This technique has been used only for the 
T] panels since the MDT are measuring in the Z-Y 
plane. On average, clusters of size two give a spa- 
tial resolution about half as for clusters of size one, 
which is below 10 mm as expected. 

6.3 TGC 

The basic structure of the TGC chambers and their 
assembly in the MS end-cap wheels is presented else- 
where |2 |. Inactive regions due to the gas-gap frame 
and the wire supports account for a loss of active area 



varying from 3% to 6% depending on the chamber 
type. In order to optimize the trigger efficiency these 
inactive regions are staggered with respect to the tra- 
jectory of high momentum muons produced at the 
IR In the active area the TGC wires are expected to 
have an efficiency of more than 98%. For the cos- 
mic ray run 91060 the trigger logic required a coin- 
cidence of 3 out of 4 layers in the doublet chambers 
(referred to as TGC2 and TGC3 as in Figure (TJ. In 
evaluating the detector efficiency one has to take into 
account the trigger bias and the fact that cosmic rays 
are non-pointing to the IP, asynchronous, and do not 
only consist of single muons but also of extended 
showers. 

To evaluate the efficiency of a layer in the dou- 
blet chambers, it is required that there is one and only 
one hit in each of the other three layers and that these 
three hits are associated to the current BC. This is in- 
tended to remove high multiplicity events (showers) 
and out-of-time tracks. As a result of this selection, 
the 3 out of 4 trigger condition is satisfied indepen- 
dently of the presence of a hit in the layer under eval- 
uation. The efficiency of this layer is thus determined 
in an unbiased way. 

A similar procedure is used for the triplet cham- 
bers (TGCl). When evaluating, the efficiency of a 
layer, it is required i) that the other two layers sat- 
isfy the 2 out of 3 trigger coincidence and ii) that the 
line joining the two hits (track) crosses the layer in 
its active area. 

In both cases, the layer under test is considered 
efficient if there is at least one hit associated to any 
of the previous, current or next BC. Figure|2T|on the 
left shows an efficiency map in the wire-strip (J]-(^) 
plane, and on the right its rj projection, i.e. the strip 
efficiency. Some inactive regions are clearly visible: 
the bands in Figure |2T}Left indicate the location of 
the wire supports. 

The overall efficiency, including the inactive re- 
gions, is evaluated for a fraction of TGC layers 
(about 40% of TGC doublet layers) by requiring that 
a muon track crosses the layer under test at least 10 
cm away from its edge. The muon track is defined 
using MDT hits combined with TGC hits in the lay- 
ers that are not under evaluation. Figure|22]shows the 
distribution of the wire efficiency for different values 
of high voltage setting: 2650, 2750, 2800 and 2850 
V. The average value of the efficiency, at the nominal 
voltage of 2800 V, is 92% consistent with the local 
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Fig. 16. Single tube tracking efficiencies witii a 5(7 association cut, as explained in the text, for a BML chamber. The 
plot on the right shows an expanded view in the region where two tubes with the wire disconnected were found. 
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Fig. 17. Distribution of the tracking efficiency, with a 5cT 
hit association cut, for ~81K drift tubes in the barrel MDT. 
About 0.2% of tubes were not working and have efficiency 
compatible with zero. 



Fig. 18. Distribution of the time difference between the 
two RFC layers of a pivot plane in the ^ projection. The 
binning of the histogram corresponds to the strip read-out 
time encoding of 1/8 of BC. 



7 MDT optical alignment 



efficiency measured as explained above and the con- 
tribution from inactive-regions. 



The design transverse momentum resolution at 1 TeV 
of the MS is about 10%, this translates into a sagitta 
resolution of 50 [xm. The intrinsic resolution of MDT 
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Fig. 19. Distribution of the average efficiency for RPC 
of the Middle stations for run 91060. The two distribu- 
tions refer to two different triggers: RPC trigger (full line, 
91.33% peak efficiency) and calorimeter trigger (dashed 
line, 92.0% peak efficiency). Both distributions are nor- 
malized to unit area. The measured efficiency is lower 
than expected mainly because the read-out timing was 
still not optimal. 



Fig. 20. Distribution of the spatial resolution provided by 
the rj strips for RPC of the Middle stations. The spatial 
resolution is divided by the strip pitch. The distributions 
are normalized to unit area. 
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Fig. 21. Left: efficiency map for a TGC chamber layer. The horizontal axis is the strip channel and the vertical axis is the 
wire channel. Right: efficiency projection to the strip channels. Observed efficiency drops are consistent with the wire 
support locations. 



chambers contributes a 40 |a.m uncertainty to the track 
sagitta, hence other systematic uncertainties (align- 
ment and calibration) should be kept at the level of 
30 |xm or smaller. Since long-term mechanical sta- 
bility in a large structure such as the MS cannot be 



guaranteed at this level, a continuously running align- 
ment monitoring system |21 1 has been installed. This 
system is based on optical and temperature sensors 
and detects slow chamber displacements, occurring 
at a timescale of hours or more. The information from 
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Fig. 22. Distribution of the TGC wire efficiency of individual layers for different high voltage values. The distribution 
for 2800V, the nominal voltage in 2008, was obtained with run 91060. 



the alignment system is used in the offline track re- 
construction to correct for the chamber misalignment. 
No mechanical adjustments were made to the cham- 
bers after the initial positioning. The system con- 
sists of a variety of optical sensors, all sharing the 
same design principle: a source of light is imaged 
through a lens onto an electronic image sensor acting 
as a screen. In addition to optical position measure- 
ments, it is also necessary to determine the thermal 
expansion of the chambers. In total, there are about 
12000 optical sensors and a similar number of tem- 
perature sensors in the system. Optical and temper- 
ature sensors were calibrated before the installation 
such that they can be used to make an absolute mea- 
surement of the chamber positions in space, rather 
than only following their movements with time rela- 
tive to some initial positions. 



7.1 End-cap chamber alignment 

The end-cap chambers and their alignment sys- 
tem 1 22 1 were installed and commissioned during 
2005-2008, and continuous alignment data-taking 
with the complete system started in summer 2008. 
After commissioning, more than 99% of all align- 
ment sensors were operational, and only a small 
number failed during the data-taking in 2008. The 
effect of the missing sensors on the final alignment 
quality is negligible. 

The position coordinates, rotation angles, and de- 
formation parameters of the chambers are determined 
by a global minimization procedure. The total x^, 
as well as the contributions of the individual sensor 
measurements to the (pulls) can be used to es- 
timate the alignment quality from the internal con- 
sistency of the fit. If the observed sensor resolutions 
agree with the design values, one expects approxi- 
mately x^/^df — 1 and a pull distribution with zero 
mean and unit RMS width. Figure |23] shows the ob- 
served and expected pull distributions in the end- 



41 



»,2000 
S18OO 
1600 
1400 
1200 
1000 
800 
600 
400 
200 



A) 



Mean 


0.004 


RMS 


1.014 


Constant 


1275. 


Mean 


0.002 


Sigma 


0.904 




ATLAS 



-2-101 2 
Endcaps A-t-C - all sensors 

Mean 
RMS 
Constant 
Mean 
Sigma 



pull 



0.000 
0.846 
1452. 
0.000 
0.836 




ATLAS 



-2-101 2 
simulation - all sensors 



Fig. 23. The observed (A, from data) and expected (B, 
from simulation) pull distributions for the end-caps, as- 
suming design resolution for all sensor types. Correlations 
and weakly constrained degrees of freedom cause the ex- 
pected pull distribution to have RMS width below unity. 
The observed ;|^^/ndf from the fit on data is 1.4, while the 
one from simulation is 1.0. 



caps, obtained by assuming the design resolutions 
for all sensor types. 

In a second step, the assumed sensor resolutions 
are adjusted until the observed pull distributions, bro- 
ken down by sensor type, agree with the expected 



distribution. This yields the observed sensor resolu- 
tions, which are used as input to a Monte Carlo sim- 
ulation of the alignment system. The simulation pre- 
dicts a sagitta accuracy due to alignment of about 
45 [a.m, which is close to the design performance. 

Validating the alignment as reconstructed from 
the optical sensor measurements requires an external 
reference. During chamber installation, surveys of 
the completed end-cap wheels were done using pho- 
togrammetry, and the chamber positions measured 
with the alignment system agreed with the survey 
results within 0.5 mm, the quoted accuracy of the 
survey. While establishing confidence in the optical 
system, the full validation of the alignment can only 
be done with tracks. Thus, cosmic muons recorded 
during magnet-off running were used to cross-check 
the alignment provided by the optical system. 

For a perfect alignment, the reconstructed sagitta 
of straight tracks should be zero for each EI-EM- 
EO measurement tower (note that, when averaged 
over many towers, the mean value can be acciden- 
tally compatible with zero despite single towers be- 
ing significantly misaligned). For cosmic muons, the 
observed width of the sagitta distribution is domi- 
nated by multiple scattering. A shifted and/or broad- 
ened distribution would indicate imperfections of the 
alignment. Triplets of track segments were selected 
in the EI-EM-EO chambers, requiring the three seg- 
ments to be in the same sector and assigned to the 
same reconstructed track. Some segment quality cuts 
were applied for this analysis: i) ;if^/ndf < 10 and 
at most one expected hit missing per chamber; ii) 
the angle between the segments and the straight line 
joining the segments in EI and EO was required to 
be smaller than 5 (50) mrad in the precision (second) 
coordinate; Hi) at least one trigger hit in the second 
coordinate was required to be associated to the track. 
A total of 1700 segment triplets passing the cuts were 
selected in run 91060. 

Figure |24] left shows, for the two end-caps, the 
observed sagitta distribution before and after apply- 
ing alignment coiTections (i.e. the chamber posi- 
tions, rotations, and deformations as determined by 
the optical system, as well as a correction for the 
gravitational sag of the MDT wires). Figure l24l right 
shows the corresponding difference in angle in the 
precision coordinate between each of the segments 
and the track (the straight line joining the EI and EO 
segments). For the distribution on the right, the cut at 



42 



225 
200 
175 
150 
125 
100 
75 
50 
25 



A) 



^ nominal geometry 

□ alignment corrections 



Constant 
Mean 
Sigma 1 
Sigma 2 



-0.033 
0.928 
3.846 



ATLAS 



-25 -20 -15 -10 



B) 



S nominal geometry 
□ alignment corrections 



10 



Constant 
Mean 
Sigma 1 
Sigma 2 



15 20 25 
sagitta [mm] 

1 01 .6 / 47 
284.6 
-0.033 
0.589 
2.297 



ATLAS 



10 -8 -6 -4 -2 2 4 6 8 10 

segment-track angular difference [mrad] 



Fig. 24. A): measured sagitta distribution for the two end- 
caps. Thie cross-iiatciied iiistogram shiows tlie sagitta be- 
fore alignment corrections, thus reflecting the accuracy 
of chamber positioning. The filled histogram shows the 
sagitta after applying alignment corrections, the curve is 
the fit of a double-Gaussian function, each Gaussian con- 
taining 50% of the events. B): measured angle in the pre- 
cision coordinate between the segments and the track to 
which they are associated. 



5 mrad was omitted. The improvement in both vari- 
ables is clearly visible, the mean value of the cor- 
rected sagitta distribution as obtained from the fit 
with a double-Gaussian function is (—33 ±42) |j.m 
and thus perfectly compatible with zero within the 
45 [a.m error estimated above from the internal con- 
sistency of the alignment fit. The width of the cor- 



rected sagitta distribution agrees approximately with 
expectations for the typical energies of triggered cos- 
mic muons. The width of the corrected angle distri- 
bution, on the other hand, is about twice as large as 
expected. This is mainly a consequence of the addi- 
tional time jitter of MDT measurements described in 
Section |5] which deteriorates the segment resolution. 

For the two end-caps separately, the mean value 
of the sagitta distribution is (— 30±61) \xm in side A 
and (—37 ± 57) [a.m in side C. The sign of the sagitta 
is defined in such a way that most of the conceivable 
systematic errors would cause deviations from zero 
with the same sign in side A and side C. The analysis 
is limited by statistics even though it uses a signifi- 
cant fraction of the full 2008 data sample. Breaking 
it down further to the level of sectors, or even to pro- 
jective towers (where the best sensitivity is obtained) 
would require significantly more data. 

The cross-check with straight tracks confirms that, 
with the limitations of the analysis, the chamber po- 
sitions given by the optical alignment system are within 
the estimated sagitta uncertainties, indicating that the 
optical system works properly. The design accuracy 
has nearly been reached in the end-caps. It also shows 
that the system produces a reliable estimate of the 
uncertainty of the alignment corrections. 



7.2 Barrel chamber alignment 

The installation and commissioning of the barrel op- 
tical system |23| began in 2005 and continued to- 
gether with the installation of the chambers until 2008. 
At the time of recording cosmic ray data, the bar- 
rel optical system was fully installed and 99.7% of 
the sensors were functioning correctly. Table |4] sum- 
marizes the status of the 5800 installed sensors. The 
complete system is read out continuously, at a rate of 
one cycle every 20 minutes. The readout was func- 
tioning correctly during the complete period of ac- 
quisition of cosmic ray data. 

The alignment reconstruction consists in deter- 
mining the chamber positions and orientations (re- 
ferred to as "alignment corrections") from the opti- 
cal sensor measurements. This requires the precise 
knowledge of the positions of the sensors with re- 
spect to the MDT wires. To this purpose, the optical 
sensors were calibrated before installation and their 
mechanical supports were glued with precise tools 
onto the MDT tubes. However, the original design of 
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Type 


Total 


Working 


Broken 


Projective 


117 


117 





Axial 


1036 


1031 


5 


Praxial 


2010 


2008 


2 


Reference 


256 


253 


3 


CCC 


260 


260 





BIR-BIM 


32 


32 





Inplane 


2110 


2101 


9 


Total 


5817 


5798 


19 


% 




99.7 


0.3 



Table 4. Status of the barrel optical system in Fall 2008. 
No data were recorded during this period from the "bro- 
ken" sensors. Naming and functions of the different sen- 
sors are detailed in reference 1231 . 



the baiTel optical system suffered from a few eiTors 
that eventually degraded the precision of the align- 
ment coiTections. Furthermore, the only devices giv- 
ing projective information in the Small sectors are 
the CCC sensors which are designed to provide 1 
mm accuracy. The alignment of the chambers of the 
Small sectors is, by design, based on tracks that cross 
the overlap region between the Small and the Large 
sectors. However, the statistics obtained in cosmic 
runs was not sufficient to perform a precise check of 
this method. 

The alignment coiTections discussed here cover 
the nine upper sectors (1 to 9). The complete period 
of cosmic data taking was divided in intervals of 6 
hours, and alignment coiTections were reconstructed 
using the sensor measurements recorded in that inter- 
val. This provided data for monitoring of significant 
movements of the MS, e.g. when the magnetic field 
in the toroids was switched on. 

The barrel alignment reconstruction is based on 
the minimization of a x^, whose inputs are, for each 
optical sensor /: 

- the recorded response r,-; 

- a model m,- (a) , representing the predicted response 
of sensor / with respect to the alignment correc- 
tions a; 

- the error (7, , the estimated uncertainty of the model 
m, 

The critical part is the model m,-, as it combines all 
the knowledge of the precise geometry of the opti- 
cal sensors and their calibration. The free parame- 
ters in the fit are the alignment corrections a, and in 
some cases additional parameters used to model the 



effect of imprecise sensor positioning or of an incor- 
rect calibration. For all these additional parameters 
appropriate constraints are included in the fit reflect- 
ing the best estimates of the error contributions men- 
tioned above. Overall, 4099 parameters are fit simul- 
taneously. The total reconstruction time for the full 
barrel is less than one minute. 

Given the uncertainties introduced by the addi- 
tional parameters in the fit procedure, the strategy for 
alignment in the barrel is slightly different from the 
one in the end-cap. Dedicated runs without magnetic 
field in the toroids (but with field in the solenoid to 
tag high momentum tracks) will be used to get ini- 
tial alignment corrections with a precision of 30 [im. 
The optical alignment system is then used to moni- 
tor movements due to the switching on of the toroidal 
field and to temperature effects. The mechanical sta- 
bility of the system, in periods where the magnetic 
field was constant, is at the level of 100 |J.m, while 
movements of the magnet structures at the level of 
few mm were observed when the magnets were switched 
on and off. The optical alignment system, which con- 
tinuously monitors the position of the chambers, is 
able to follow these movements with the required ac- 
curacy. This so-called relative alignment mode has 
already been tested with success in the MS system- 
test done with a high-energy muon beam [4^1. Af- 
ter the minimization, the value obtained for ;|f ■^/ndf is 
1 .9, which shows that the sensor errors are underes- 
timated. 



7.2.1 Performance of the optical alignment in 
the barrel 

Similarly to what is done in the end-cap, an esti- 
mate of the contribution to the sagitta error due to 
the alignment system may be inferred from the 
using the following formula 



^ IdOidOj e 

where 0, are the fitted parameters and V is the 
global error matrix, of size 4099x4099, of all fit- 
ted parameters. To estimate the performance of the 
alignment system in terms of sagitta measurement, 
straight tracks, originating in the IP and crossing 
three layers of chambers, were simulated and the 
whole fit procedure was applied to these tracks. The 
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sagitta of these pseudo-tracks is a function of some 
of the ahgnment corrections, and thus the formula 
of error propagation may be used to infer the contri- 
bution of the alignment to the error of the resulting 
sagitta. This technique relies on the hypothesis that 
the errors of the optical sensors are correctly esti- 
mated, and thus that the is correctly normalized. 
As this is not the case (x^/ndf = 1.9), the results are 
only considered as a rough estimate of the optical 
alignment performance. 

The result is shown in Figure |25] The Small sec- 
tors have a significantly worse alignment than the 
Large sectors, as explained above. Conservatively, 
one can conclude that the performance of the optical 
system, in terms of sagitta precision, is ^200 [xm for 
the Large sectors, and ^1 mm for the Small sectors. 




Fig. 25. rj X (j) map of the contribution to the sagitta er- 
ror due to alignment, as estimated with the method de- 
scribed in the text. As expected from the system design, the 
Small sectors (even sector numbers) are aligned with sig- 
nificantly less precision than the Large sectors (odd sector 
numbers). 



7.2.2 Alignment with straigint tracks 

Data with the toroidal field off were used to improve 
the alignment precision in the barrel and to validate 
the alignment corrections in relative mode. The method 
is to use straight muon tracks to determine in ab- 
solute mode the initial spectrometer geometry and, 
once this geometry is determined, to use the optical 
alignment system to trace all chamber displacements 



in a relative mode. The alignment procedure with 
straight tracks is based on the so-called MILLEPEDE 
fitting method |24|. This method uses both alignment 
and track parameters inside a global fit. As a result, 
all correlations between alignment and track param- 
eters are taken into account and the alignment algo- 
rithm is unbiased. 

The track alignment algorithm has been tested 
with Monte Carlo simulations and with cosmic ray 
data. The simulation studies show that 10^ muon tracks 
with a momentum greater than 20 GeV and pointing 
to the IP are needed to align the Large sectors with a 
precision of 30 |J.m. Small sectors require five times 
more tracks than Large sectors, due to the multiple 
scattering in the toroid coils. 

Using straight cosmic muon tracks recorded in 
run 91060, a set of alignment constants has been pro- 
duced. A total of 10^ events were used corresponding 
to about 3x10^ cosmic muon tracks in each of the 
most illuminated barrel sectors. The statistical uncer- 
tainty of the sagitta using this track alignment proce- 
dure was estimated to be 30 |J.m for Large sectors. 

The data of run 91060 were processed with the 
track reconstruction software twice: /) using the op- 
tical alignment corrections and //) using the track- 
based alignment corrections. Both geometries were 
then tested by measuring the distribution of the track 
sagitta for muons crossing three chamber stations (In- 
ner, Middle and Outer). Only tracks passing close to 
the IP in the rj projection were chosen. Hits in the 
Inner and Outer chambers were fit to a straight line, 
and the distribution of the hit residuals in the Middle 
chambers was used to evaluate the sagitta. For per- 
fect alignment, the mean value of the sagitta should 
be zero for straight tracks and, to a good approxi- 
mation, the mean value of the distribution gives an 
estimate of the sagitta error. 

The results are shown in Figure|26]for the sets of 
alignment corrections; on the left for a station in a 
Large sector, on the right for a station in a Small sec- 
tor For reference, the distributions using the design 
geometry are also shown. The tails of the distribu- 
tions are due to multiple scattering of muons. In the 
Large sector station, the two distributions are almost 
identical, but the distribution with optical alignment 
is centered at ^100 [a.m. In the Small sector station, 
the distribution with the optical alignment is centered 
around 1 mm. To compare the results obtained in dif- 
ferent stations. Figure |27] shows the mean values of 
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Fig. 26. Distribution of tlie sagitta (as defined in the text) for straight tracks. A), B): using alignment corrections derived 
from the optical system only; C),D): using track-based alignment corrections. A), C): for a station in a Large barrel 
sector; B), D): for a station in a Small barrel sector, the optical system corrections of the small sectors have, by design, 
an accuracy at a level of 1 mm. In all panels, the hashed distribution is obtained using the "nominal" geometry. Mean 
and sigma in the statistical boxes refer to the distributions with alignment corrections, the peak is fitted with a gaussian 
in a ± 1 sigma interval 



the sagitta distribution for the Large upper sectors (3, 
5 and 7). One Small sector is also presented, sector 4, 
since this was illuminated with enough events during 
the same run to produce a meaningful distribution. 

The results show that the optical alignment 
system alone provides a precision at the level of 



200 |i.m. When calibrated with sufficient statistics of 
high momentum straight tracks, the optical system is 
able reach a precision of 50 [a.m. 

The sagitta resolution for runs with no magnetic 
field in the MS can be studied as a function of the 
muon momentum measured by the ID. The sagitta 
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Fig. 27. Mean value of the track sagitta distributions obtained A), B): with the optical alignment system only, C), D): 
and using the track-based alignment . A), C): for the upper Large barrel sectors. B), D):: for a Small barrel sector with 
56° <(j><79°. 



resolution as a function of the muon momentum was 
parameterized as 

where the first term Kq is due to multiple scatter- 
ing in the material of the MS, and the second term 
Ki is due to the single tube resolution and chamber- 
to-chamber alignment. These two terms have been 



already measured at the MS system test beam ID 
|5] and found to he Kq = 9 mmxGeV and Ki = 
50 |xm. A similar measurement was done with cos- 
mic muons by selecting segment triplets (Inner, Mid- 
dle and Outer station) of MS projective towers. The 
RMS of the sagitta of the Middle station segment 
with respect to the Outer-Inner straight line extrapo- 
lation has been fit in five momentum bins. The result 
is shown in Figure|28]for sector 5 (Large sector) with 
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RPC-time corrections applied in the calibration pro- 
cedure. The fitted value for the two terms is Kq ~ 
(12.2±0.7)mmxGeVand/:i = (107±21) ^m. In 
the MS the multiple scattering term is expected to be 
worse than the one measured at the test beam setup 
and larger for Small sectors due to the presence of 
the toroid coils between the Inner and Outer cham- 
bers. The value of Ki measured with cosmic muons 
in sector 5 is only about a factor two worse than that 
measured at the test beam. Several effects contribute 
to this, including alignment, chamber deformations, 
calibration and single tube resolution. Similar stud- 
ies performed for other sectors show worse results 
due to the smaller data sample available for align- 
ment and calibration. 
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Fig. 28. RMS value of the sagitta distribution in sector 5 
as a function of the muon momentum measured by the ID. 
The fit to the function described in the text is superim- 
posed. 



These preliminary studies with cosmic rays in- 
dicate that the method of track-based alignment is 
robust and with sufficient muon data from collisions 
the design alignment precision will be achieved. 

8 Pattern recognition and segment 
reconstruction 

The pattern recognition algorithm first groups hits 
close in space and time for each detector Each pat- 
tern is characterized by a position and a direction 
and contains all the associated hits. Starting from 
these patterns, the segments are reconstructed with a 



straight line fit. The Gfo-refit is applied at this stage 
and, if the Gfo-refit procedure does not converge, 
the segment parameters are computed with the tube 
to provided by the calibration with tube resolution 
increased to 2 mm. After this, a drift radius is as- 
signed to each tube with an uncertainty of 2 mm (in- 
dependent of the drift radius value) in order to keep 
high track reconstruction efficiency even in the case 
where no precise alignment constants are available. 
The minimum number of hits per segment was set to 
3 and no cuts were applied on the number of missed 
hits. 

These relaxed requirements tend to increase the 
number of fake segments while keeping a high seg- 
ment efficiency. Since cosmic ray events are quite 
clean and have low hit multiplicity this fake rate in- 
crease is not considered as a problem. On the other 
hand, a high reconstruction efficiency allows the use 
of segments to spot hardware problems in individ- 
ual chambers or in calibration or decoding software. 
Most of the fake segments are rejected at the track 
reconstruction level. 

Figure|29]shows, on the left, the number of MDT 
hits per segment for segments associated to a track. 
In the distribution clear peaks are observed at 6 and 8 
hits corresponding to the 6-layer (Middle and Outer) 
and 8 -layer (Inner) chambers. 

The efficiency of the segment reconstruction in 
run 91060 was determined in the following way. 
First, cosmic shower events are suppressed by re- 
quiring less than 20 segments in the event. Then 
a pair of segments in two MDT stations (Inner, 
Medium or Outer) are fitted to a straight line and the 
line is extrapolated to the third station. In the extrap- 
olation multiple scattering is taken into account as- 
suming a 2 GeV momentum for the cosmic muon. If 
the extrapolated line crossed the third station, a re- 
constructed segment is searched for in that station, 
but it is not required that the hits of this segment be 
associated to the muon track. The segment efficiency 
is then computed for each MDT chamber as the frac- 
tion of times a segment is found. In order to reduce 
the effect of the non-instrumented regions a fiducial 
cut in Tj was applied for both barrel and end-cap. 
Chambers that were not operational in the analyzed 
run were removed from the sample. It was not possi- 
ble to determine the efficiency for all chambers due 
to the limited coverage of the trigger for the run used 
for this analysis (fall 2008) and flux of cosmic rays. 
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Fig. 29. A): distribution of the number of MDT hits per 
segment for segments associated to a track. B): segment 
reconstruction efficiency for 322 MDT chambers. 



For tracks crossing the overlap region between two 
adjacent chambers ( Small/Large sector overlap) it 
was not required that two segments be reconstructed, 
since this may lead to a slight overestimation of the 
efficiency. 

The distribution of the segment efficiency is 
shown in Figure |29] on the right for 322 chambers 
in the barrel. The average value is 99.5% and the 
segment efficiency is uniform over the acceptance 



as shown in Table |5] In the efficiency for the bar- 
rel chambers there is a small loss of about 0.5% due 
the presence of the support structure of the ATLAS 
barrel. The Inner chambers have a slightly lower seg- 
ment efficiency due to the geometry of the trigger 
and a larger uncertainty in the track extrapolation. 
Studies on systematic effects in determining the seg- 
ment efficiency, such as its dependence on the Gfo- 
refit, on the extrapolation and on the track angle, 
show that a systematic error of ^0.5% affects the 
values of efficiency listed in Table |5] 

An alternative method to evaluate the seg- 
ment reconstruction efficiency, almost independent 
of chamber hardware problems, is described in the 
following. As in the previous case, this method can 
be used only with no magnetic field in the MS. All 
segment pairs in two different MDT stations (In- 
ner, Middle or Outer) with a polar angle difference 
smaller than 7.5 mrad are considered. The segment 
pairs are fitted to a straight line and this is extrapo- 
lated to the third MDT station. The track is kept if at 
least three hit tubes are found in the third MDT with 
a signal charge above the ADC cut and aligned with 
the track extrapolation within one tube diameter, ±3 
cm. The segment efficiency is then computed as the 
fraction of selected tracks that have a segment re- 
constructed with at least 3 hits in the identified drift 
tubes. Since the normalization already requires the 
presence of three hits in the tested MDT station, this 
segment efficiency is almost independent of local 
hardware problems. A segment reconstruction effi- 
ciency higher than 0.99 is found in all MDT stations. 

The rate of fake segments was studied with a 
random trigger. An average rate of 0.06 fake seg- 
ments per event was found with the relaxed hit as- 
sociation criteria used for cosmic muons. This rate 
is expected to be strongly reduced to about 2 x 10^^ 
if the segment reconstruction requirements are made 
to be more stringent as shown by using an alternative 
muon tracking algorithm. 



9 Track reconstruction 

The MOORE and Muonboy programs have been op- 
timized to reconstruct muon tracks originating from 
the IR To cope with the different topology of cos- 
mic ray muons they have been slightly modified as 
explained in Section[2] To mimic muons in collision 
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MDT station 


BI BM BO 


EI EM EO 


Segment efficiency 


0.987 0.992 0.996 


0.992 0.998 0.999 



Table 5. Average value of the segment reconstruction efficiency in the MDT stations. 



events, the tracks are split at their perigee (point of 
closest approach to the beam axis), giving, usually, 
two reconstructed tracks: one in the upper part of the 
MS and one in the lower part. Events with at least 
one ID track satisfying the following criteria were 
selected: 

- at least 20 hits in the Transition Radiation 
Tracker; 

- the number of hits summed over the Semicon- 
ductor Tracker (SCT) and the Pixel detector 
greater than 4 ; 

- the distance of closest approach in the transverse 
plane jc/ol and along the z axis |zo| smaller than 1 
m; ; 

- the value of the muon track ;i;^/ndf < 3 ; 

- the value of the reconstructed pseudorapidity 

l'?l<i; 

- reconstructed momentum greater than 5 GeV. 

This selection has been applied for all the stud- 
ies reported in this Section with the exception of the 
momentum resolution results. 



9.1 Resolution 
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Fig. 30. Distribution of residuals for MDT hits associated 
to a track. The residuals have been fitted with a double- 
Gaussian function with common mean. The mean value 
is 6 llm, the standard deviation of the narrow Gaussian is 
about 150 ilm and the one of the wide Gaussian is about 
700 /xm. 



The distribution of residuals for MDT hits associ- 
ated to a track is shown in Figure[30l The hit residual 
is defined as the difference between the drift radius 
measured in a tube and the distance of the track to the 
tube wire. The distribution refers only to tracks with 
MDT hits in at least three different muon stations 
(Inner, Middle and Outer) because these tracks have 
well constrained parameters and individual hits give 
a small contribution to the track parameters. The dis- 
tribution was fitted to a double Gaussian with com- 
mon mean value. The mean of the distribution was 
6 |J.m and the RMS widths was 150 |J.m for the nar- 
row Gaussian, accounting for 75% of the distribu- 
tion, and 700 |J.m for the other. When compared to 
the distribution of the segment residuals shown in 
Section|6]two additional effects contribute to the broad- 
ening of this distribution: the misalignment between 
stations and multiple scattering in the MS material. 



9.2 Efficiency 

The track reconstruction efficiency is computed as 
the fraction of events where a track is reconstructed 
in the MS top or bottom hemisphere once an ID track 
was found satisfying the selection criteria described 
above. In this case also tracks with hits in only two 
out of three MDT stations (Inner, Middle or Outer) 
are accepted, even if these tracks have a worse mo- 
mentum resolution than tracks reconstructed in three 
stations. About 15% of the selected tracks are in this 
category. In addition, to compute the track efficiency 
in the top (bottom) hemisphere, a momentum cut of 5 
GeV (9 GeV) on the ID track is applied. The result is 
shown in Figure|3T|as a function of the pseudorapid- 
ity of the ID track, for the top and bottom hemisphere 
separately. 
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Fig. 31. Track reconstruction efficiency as a function of 
pseudorapidity. The loss in efficiency in the region near 
1 77 1 = is due to the loss of acceptance for detector ser- 
vices. The presence of a track measured in the ID with 
1 77 1 < 1 is required. 

The value of the efficiency, integrated over the 
Tj acceptance, is 94.9% for the top and 93.7% for 
the bottom hemisphere respectively. If the four cen- 
tral bins are removed in Figure [3T| the efficiency in- 
creases to 98.3% and 96.3% respectively. The statis- 
tical error on these values is below 0.1%. The lower 
efficiency in the central detector region, around | rj | = 
0, is due to the presence of the main ATLAS service 
gap while lower efficiency in the Bottom hemisphere 
is explained by the uninstrumented regions occupied 
by the support structure of the ATLAS barrel. 



9.3 Momentum measurement 
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Fig. 32. A): Distribution of momentum of cosmic muons 
as measured at the MS entrance for the upper and lower 
hemispheres. The difference between the two distributions 
is due to the ID track momentum cut of 5 GeV. B): same 
distributions with track momentum extrapolated to the IP. 



The momentum of cosmic muons was measured in 
runs with magnetic field. The momentum measure- 
ment can be defined at the MS entrance or at the 
point of closest approach to the IP. In the second 
case, for tracks crossing the ID, a correction was 
made for the energy loss in the calorimeters. This 
correction is based on the average energy loss com- 
puted by the track extrapolation algorithm and is on 



average 3.1 GeV for muons pointing to the interac- 
tion region with a distance of closest approach of 
|t/o| < 1 m and |zo| < 2 m. 

The distribution of momentum at the MS entrance 
is shown in Figure [321-left for the top and bottom 
hemispheres separately. The difference between the 
two distributions is due to the ID track momentum 
cut of 5 GeV that translates in a different momen- 
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turn cut-off in the two MS hemispheres, since cosmic 
muons are directed downwards. The same distribu- 
tion extrapolated to the perigee is shown on the right 
side of Figure [32l demonstrating that the correction 
for the energy loss in the calorimeter removes the 
offset. 
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Fig. 33. A): Number of MDT hits on track. B): momentum 
difference between momenta measured by the MS in the 
top and bottom hemispheres for cosmic muons. The mo- 
menta are expressed at the MS entrance and only tracks 
with momenta bigger than 15 GeV are considered. The 
mean value of 6.3 GeV is due to the energy loss in the 
calorimeter material. 



The distribution of the number of MDT hits as- 
sociated with a track is shown on the left side of Fig- 
ure |33] For this plot tracks measured in three MDT 
stations have been selected. A clear peak around 20 
hits is visible (8 tubes in the Inner stations, 6 in the 
Middle and Outer stations). 

In events with tracks that cross the whole MS, 
the track is split at the perigee and the two indepen- 
dent momentum measurements, in the top and bot- 
tom hemisphere, can be compared. Figure [33] (right 
side) shows the distribution of the difference of the 
two momentum values, top-bottom measured at the 
MS entrance, for tracks with momenta greater than 
15 GeV. In this case the muons cross the calorimeter 
twice and the energy loss is twice the value quoted 
above, in good agreement with the 6.3 GeV mean 
value of the distribution. 

The MS momentum resolution has been estimated 
by comparing for each cosmic muon the two inde- 
pendent measurements in the top and bottom hemi- 
spheres. In order to increase the available statistics 
no requirements on the presence of ID tracks were 
applied in this study. Only events with at least two 
reconstructed tracks in the MS are considered. Each 
track is required to have: 

- at least 17 MDT hits, of which at least 7 in the 
Inner and 5 in the Middle and Outer stations of 
the same sector; 

- at least 2 different layers of RPC with a hit in the 
projection; 

- polar angle 65° < < 115°; 

- distance of closest approach to the IP |t/o| < 1 m 
and |zo| < 2 m ; 

- polar and azimuthal angles of the MS track pair 
agree within 10 °. 

About 19K top-bottom track pairs were selected 
in this way. For each track the value of transverse 
momentum was evaluated at the IP. The difference 
between the two values divided by their average 

^PT _ 2 PTup — PTdown 



PT 



PTiip + PTdo 



was measured in eleven bins of pr- Since the cosmic 
muon momentum distribution is a steep function (see 
Figure |32]l. the pj value of each bin was taken as the 
mean value of the distribution in that bin. 

The distribution of Apr/pr was fitted in each 
bin with a double-Gaussian function with common 
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mean value. The narrow Gaussian was convoluted 
with a Landau function to account for the distribu- 
tion of energy loss in the calorimeter For pr < 10 
GeV the normalizations of the two Gaussians were 
constrained such that 95% of the events are in the 
narrow Gaussian. Above 10 GeV this constraint was 
lowered to 70%. The mean value is representative 
of the difference in the transverse momentum scale 
between the two MS hemispheres. The RMS of the 
narrow Gaussian plus the width of the Landau, di- 
vided by a/2, is taken as an estimate of the transverse 
momentum resolution for each pj bin. The Landau 
width is added linearly to the narrow Gaussian RMS 
since the two quantities are strongly correlated. 

The distribution of Apr/pr is shown in Figure[34 
for all pt bins together with the fitted function. For 
the eleven bins the fit probability is in the range be- 
tween 45% and 99%, showing that the chosen parametriza 



0.6) X 10^^ GeV^'. The second uncertainty, due the 
systematics of the bin-by-bin fit method, was evalu- 
ated by changing the fitting assumptions as explained 
above. The expected values for these parameters were 
computed in reference |3| on the basis of an analytic 
calculation of the pj resolution that takes into ac- 
count the detailed description of the material in the 
MS, the single tube resolution, the alignment accu- 
racy and the magnetic field map. The values obtained 
for the barrel MS were: Pq = 0.35 GeV, Pi = 0.035 
and P2 = 1.2 x 10""* GeV^^ The result is in fair 
agreement with the expected values for the first two 
terms, while the intrinsic term is worse. The differ- 
ence has been investigated to trace the effects that 
contribute to worsen the resolution as determined with 
cosmic muons. 



tion is a good representation of the data distribution. 

Different fits have been done to study the system- 
atics of the mean and RMS value, i) The constraint 
between the two Gaussian areas has been changed 
by ±10%. ii) A double Gaussian with common mean 
and asymmetric fit range, with the fit range reduced 
to two standard deviations on the positive side to 
avoid the energy loss tail. Hi) A fit with two indepen- 
dent Gaussians with no range constraint. The result 
is that the estimated resolution is quite independent 
of the fit assumptions. The variation of the fit reso- 
lution ranges between 0.5% at low pr up to a maxi- 
mum of 1 % in the highest momentum bin. 

The fit mean values indicate that the pr scales 
in the two MS hemispheres are in agreement within 
1%, or better The relative pr resolution, Opj/pr = 

o{ApT / Pt) / is shown in Figure[35] for the two 
main muon reconstruction algorithms ||7][8l, as a func- 
tion of the transverse momentum. The two results are 
consistent taking into account the independent statis- 
tical uncertainties. 

The resolution function can be fitted with the sum 
in quadrature of three terms, the uncertainty on the 
energy loss corrections Pq, the multiple scattering 
term Pi, and the intrinsic resolution P2. 



First, more than 70% of the track pairs consid- 



PT 



PT 



Pi X PT 



The result of the fit is shown in Figure |35] The val- 
ues of the parameters are: Pq = 0.29 ± 0.03 ± 0.01 
GeV, Pi = 0.043 ± 0.002 ± 0.002, P2 = (4. 1 ± 0.4 ± 



ered in the analysis are in the Large sectors 5-13. 
At high Pt the momentum resolution in the barrel 
Large sectors is worse than in Small sectors because 
the field integral is smaller (see Figure [T]). Instead, 
in the low pT region dominated by multiple scatter- 
ing the resolution in Large sectors is better since the 
magnet coils are in the Small sectors. 

Second, the single tube resolution is affected by 
imperfect calibrations and the additional time jitter is 
not completely recovered by the Gfo-refit (see Fig- 
ure[T2]|. Part of the tracks in the sample contain seg- 
ments with a badly converging Gfo-refit. As a cross- 
check, all the tracks with bad convergence were re- 
moved and the analysis was repeated. The intrinsic 
term decreased by about 30%. 

Third, the alignment in many sectors of the MS 
is still not at the required level due to the limited 
statistics of straight tracks in the cosmic ray data 
sample. Last, several other effects that contribute to 
resolution have not been removed, such as chamber 
deformations (due to temperature effects), wire sag- 
ging (particularly important in large chambers), sin- 
gle chamber geometrical defects. Each of these ef- 
fects contribute to worsening the resolution and can 
be removed with dedicated software tools. At the 
present stage of commissioning, the momentum res- 
olution is close to the design value for pT < 50 GeV, 
but is not as good for higher momenta. 
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Fig. 34. Distributions of Apj/pj in the eleven pj bins. Fits to the function described in the text are superimposed. 
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Fig. 35. Transverse momentum resolution evaluated with the top-bottom method explained in the text as a function of 
pj, barrel region only (|t)| < 1.1). The fit to the three resolution parameters as described in the text is superimposed. 



10 Summary 

The data collected in several months during the 2008- 
2009 cosmic ray runs have been analyzed to assess 
the performance of the Muon Spectrometer after its 
installation in the ATLAS experiment. Parts of the 
detector, the Small Wheels in front of the end-cap 
toroids, were installed during the runs and the com- 
missioning of the many detectors was proceeding while 
debugging the data acquisition and the data control 
systems. The detector coverage during most of the 
run period was higher than 99%, with the excep- 
tion of the RPC chambers which were still under 
commissioning. For this detector subsystem the cov- 
erage steadily improved during the commissioning 
runs reaching more than 95% in Spring 2009. Re- 
sults on several aspects of the Muon Spectrometer 
performance have been presented. These include de- 
tector coverage, efficiency, resolution and relative tim- 
ing of trigger and precision tracking chambers, track 
reconstruction, calibration, alignment and data qual- 
ity. 

Finally, with data collected when the magnetic 
field was on, a first estimate of the spectrometer mo- 



mentum resolution was obtained. Efficiency and res- 
olution of single elements have been measured for 
MDT, RPC and TGC chambers and were found in 
agreement with results obtained previously with high- 
momentum muon beams. The trigger chamber tim- 
ing has been adjusted with enough precision to guar- 
antee that the interaction bunch crossing can be iden- 
tified with a minimal number of failures. The muon 
trigger logic, based on fast tracking of pointing muons 
has been extensively tested in the regions of the de- 
tector with good cosmic ray illumination. A slight 
deterioration of the MDT spatial resolution, com- 
pared to test beam results, was observed, which can 
be understood in terms of an additional time jitter 
due to the asynchronous timing of cosmic muons and 
to their non-pointing geometry. These effects were 
partially removed, modifying the track reconstruc- 
tion programs with dedicated algorithms. Allowing 
for an increase of the single hit resolution, to cope 
with these effects, the track segment efficiency in in- 
dividual chambers was found to be satisfactory and 
uniform over the large number of chambers. 

The performance of the end-cap and baiTel op- 
tical alignment systems have been measured using 
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cosmic muon tracks with no magnetic field. The re- 
sults demonstrate that the end-cap optical system is 
able to provide the required precision for chamber 
alignment. The design of the alignment system in 
the barrel requires additional constraints provided by 
straight tracks. The method has been tested with good 



in the collaborating Institutions without whom the 
ATLAS detector could not have been built. Further- 
more we are grateful to all the funding agencies which 
supported generously the construction and the com- 
missioning of the ATLAS detector and also provided 
the computing infrastructure. 



results, but is limited by the statistics of high-momentum -^1^ ^^LAS detector design and construction has 



muons with the required pointing geometry. 

With the geometry corrections provided by the 
alignment system, tracks in projective geometry 
were reconstructed in the barrel showing that the re- 
construction efficiency is uniform over the entire ac- 
ceptance and that the sagitta error is in agreement 
with the detector resolution, the alignment precision 
and the effect of multiple coulomb scattering. 

Finally with magnetic field, tracks crossing the 
whole spectrometer were used to obtain two inde- 
pendent measurements of the momentum. The mo- 
mentum resolution was evaluated using the two val- 
ues in the top and bottom part of the detector and 
the results were analyzed, fitting the distribution of 
the difference as function of the momentum. Taking 
into account the momentum spectrum, the multiple 
scattering in the spectrometer and the energy loss in 
traversing the calorimeters, the momentum resolu- 
tion is in good agreement with results from simu- 
lation for transverse momenta smaller than 50 GeV. 
The statistics of high-momentum pointing tracks lim- 
its the accuracy of the individual chamber calibration 
and the precision of the alignment. At higher mo- 
menta, these limitations result in degraded momen- 
tum resolution. 

During the long period of commissioning with 
cosmic rays it was possible to optimise the perfor- 
mance of the various hardware and software elements 
and to reach a level of understanding, such that we 
can consider the Muon Spectrometer to be ready to 
efficiently detect muons produced in high-energy proton' 
proton colUsions 



taken about fifteen years, and our thoughts are with 
all our colleagues who sadly could not see its final 
realisation. 
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